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I.  Introduction 


The  work  performed  in  this  effort  is  directed  towards  the  design, 

analysis,  and  prediction  of  electromagnetic  compatibility  (EMC)  in 

3 

command,  control,  and  communications  (C  )  systems.  This  was  accomplished 
by  closing  some  basic  technological  gaps  in  two  important  areas.  The 
first  area  dealt  with  electromagnetic  (EM)  effects  on  microelectronic  circuits 
while  the  second  was  concerned  with  EM  coupling  and  testing. 

To  obtain  insight  into  the  response  of  microelectronic  circuits 
to  electromagnetic  interference  (EMI),  EMI  effects  were  studied  in 
several  types  of  operational  amplifiers  and  in  a  differential  line  receiver. 
Laboratory  measurements  were  performed  on  the  operational  amplifiers 
whereas  a  computer  simulation  was  made  of  the  differential  line  receiver. 

Attention  was  also  devoted  to  the  problem  of  analyzing  and  predic¬ 
ting  EMI  in  microelectronic  circuits.  It  was  shown  that  the  task  of 
modeling  complicated  digital  circuits  for  EMI  effects  can  be  greatly 
simplified  through  the  use  of  macromodels.  In  particular,  macromodels 
were  developed  for  operational  amplifiers,  a  NAND  gate,  and  a  four  bit 
comparator.  In  addition,  the  probabilistic  approach  to  EMI  was  inves¬ 
tigated.  This  included  an  experimental  study  of  the  statistical  varia¬ 
tions  of  EMI  in  operational  amplifiers,  a  probabilistic  analysis  of  EMI 
induced  propagation  delays  in  combinational  digital  circuits,  and  the 
analysis  of  EMI  in  an  experimental  system  containing  both  a  random 
coupling  mechanism  and  a  random  second-order  rectification  effect. 

To  eliminate  gaps  in  the  area  of  EM  coupling  and  testing,  studies 
were  carried  out  in  1)  coupling  of  EM  fields  into  transmission  lines, 


2)  crosstalk  between  shielded  cables,  3)  crosstalk  in  twisted  wire 


circuits.  4)  crosstalk  in  balanced  twisted  pairs,  and  5)  modeling  of 
transmission  lines.  Finally,  recommendations  were  made  for  revisions 
to  the  wire-to-wire  coupling  subroutine  in  IEMCAP. 

Results  of  these  efforts  are  reported  in  six  volumes  whose  con¬ 
tents  are  briefly  summarized  below. 

Volume  IA:  "Nonlinear  Macromodel  of  the  Bipolar  Integrated  Circuit 
Operational  Amplifier  for  Electromagnetic  Interference 
Analysis"  by  Gordon  Chen  and  Dr.  James  Whalen  [1], 

This  report  describes  a  small-signal  nonlinear  maciomodel  for  the 
bipolar  IC  Op  Amp  which  can  be  used  successfully  for  RFI  analysis.  Two 
BJT's  are  adequate  to  characterize  the  complete  IC.  The  effectiveness 
of  the  small-signal  nonlinear  macromodel  is  verified  by  the  agreement 
between  the  NCAP  computer  program  predictions  and  experimental  results 
for  the  second  order  transfer  function  of  the  741  Op  Amp  unity  gain  buffer 
amplifier.  Good  agreement  is  also  obtained  between  NCAP  calculations 
based  on  the  macromodel  and  the  full  model.  By  using  the  Op  Amp  macro¬ 
model,  a  saving  of  nearly  an  order  of  magnitude  in  computer  costs  is 
achieved.  Not  only  does  the  macromodel  conserve  computer  resources,  it 
makes  it  possible  for  the  EMC  engineer  to  analyze  electronic  systems 
involving  many  Op  Amps. 

Another  important  aspect  of  the  small-signal  nonlinear  macromodel 
is  that  it  is  a  global  model  which  can  be  used  for  all  Op  Amp  types. 


2 


The  two  macromodel  transistors  should  be  the  same  type  as  those  at 
the  input  of  the  IC  Op  Amp  being  modeled.  The  need  for  laboratory 
characterization  of  the  individual  transistors  in  each  new  bipolar 
Op  Amp  appear  not  to  be  required.  A  combination  of  manufacturer's 
specifications  and  previous  IC  BJT  results  can  be  used.  This  was 
demonstrated  by  the  successful  application  of  the  nonlinear  macromodel 
to  the  new  LM.IO  IC  Op  Amp  for  RFI  analysis  by  NCAP,  without  resorting 
to  full  model  analysis,  nor  experimental  device  characterization. 

However,  the  accuracy  of  the  nonlinear  macromodel  relics  on  the 
inclusion  of  capacitors  to  represent  the  effects  of  substrate  para¬ 
sitic  capacitances.  The  proper  selection  of  substrate  capacitors  does 
not  necessarily  follow  the  structural  IC  parasitic  capacitances.  The 
effects  caused  by  these  capacitances  are  simulated  by  four  capacitors 
in  the  macromodel.  The  best  values  for  these  capacitors  were  found  by 
comparing  NCAP  predicted  values  with  experimental  values  for  a  second 
order  transfer  function  of  the  Op  Amp  in  an  unity  gain  buffer  amplifier. 
Thus  a  parametric  fitting  is  involved  when  a  new  Op  Amp  is  modeled. 

That  the  combined  nonlinearities  of  tens  of  BJTs  in  a  linear  IC 
Op  Amp  can  be  represented  quite  well  by  a  pair  of  BJTs  in  a  macromodel 
is  explained  by  the  cascading  theory  of  nonlinear  transfer  functions. 

The  essential  point  appears  to  be  that  the  second  stage  (gain-stage) 
contribution  to  the  Op  Amp  second  order  transfer  function  is  reduced 
greatly  by  the  internal  feedback  capacitor  of  the  second  stage.  Also 
it  was  argued  that  the  input  stage  which  consists  of  two  cascode  stages 


in  a  differential  pair  arrangement  can  be  represented  by  two  common  emitter 
stages  in  a  differential  pair  configuration  used  in  the  Op  Amp  macromodel. 
The  validity  of  the  nonlinear  macromodel  is  therefore  supported  both  by 
experimental  verification  and  by  network  theory. 

Volume  IB:  "Radio  Frequency  Interference  Demodulation  Effects  in 

Operational  Amplifier  Circuits"  by  Yue-Hong  Sutu  and  Dr. 

James  Whalen  [2]. 

Since  monolithic  integrated  operational  amplifiers  (op  amps)  are 
important  building  blocks  in  today's  control  and  communication  systems, 
a  series  of  investigations  are  carried  out  to  determine  RFI  effects  in 
analog  circuits  using  op  amps  as  active  devices.  The  specific  RFI  effect 
investigated  is  how  amplitude-modulated  (AM)  RF  signals  are  demodulated 
in  op  amp  circuits  to  produce  undesired  low  frequency  responses  at  the 
AM-modulation  frequency.  Subsequently,  the  undesired  demodulation  re¬ 
sponses  may  be  processed  in  a  manner  similar  to  a  desired  low  frequency 
signal  by  the  low  frequency  components  that  follow  the  op  amp.  The  un¬ 
desired  demodulation  responses  are  shown  to  be  characterized  by  second- 
order  nonlinear  transfer  functions. 

Four  representative  op  amp  types  investigated  are  the  741  bipolar 
op  amp  with  its  conventional  npn  input  transistors,  the  LM10  bipolar 
op  amp  with  its  less  conventional  pnp  input  transistors,  the  LF355 
JFET-Bipolar  op  amp  with  its  JFET  input  transistors,  and  the  CA081  MOS- 
Bipolar  op  amp  with  its  MOSFET  input  transistors. 


Two  op  amp  circuits  are  investigated.  The  first  op  amp  circuit 
configuration  is  the  unity  voltage  gain  buffer  circuit  known  as  a 
voltage  follower.  The  circuit  configuration  is  also  called  noninverting 
because  both  the  intended  signal  and  the  RFI  signal  are  injected  into 
the  noninverting  input.  The  second  circuit  is  an  inverting  op  amp  con¬ 
figuration.  The  intended  signal  voltage  gain  is  determined  by  the  ratio 
of  the  feedback  resistor  R2  and  the  input  resistor  Rl.  The  investigation 
includes  several  combinations  of  Rl  and  R2  and  also  the  effects  of  an 
RFI  suppression  capacitor  in  the  feedback  path. 

Approximately  30  units  of  each  op  amp  type  are  tested  to  determine 
the  statistical  variations  of  RFI  demodulation  effects  in  the  two  op 
amp  circuits.  Statistical  parameters  such  as  means  and  standard  devia¬ 
tions  provide  quantitative  measures  on  variations  in  op  amp  RFI  suscepti¬ 
bility. 

The  Nonlinear  Circuit  Analysis  Program  NCAP  is  used  to  simulate 
the  demodulation  RFI  response,  in  the  circuit  simulation,  the  op  amp 
is  replaced  by  its  incremental  macromodel.  Values  of  macromodel  para¬ 
meters  are  obtained  from  previous  investigations  and  from  manufacturer's 
data  sheets.  Sensitivity  analysis  of  the  effects  of  variations  of 
model  parameters  upon  RFI  susceptibility  is  performed  for  the  unity 
gain  buffer  circuit  in  which  critical  parameters  were  identified.  The 
effects  of  the  RFI  suppression  capacitor  in  the  inverting  op  amp  circuit 
is  also  simulated.  The  simulated  results  are  compared  to  measurements 
and  satisfactory  agreement  is  achieved. 


Volume  IIA:  "Macromodels  of  Op  Amps  for  CADA  Applications"  by 


Dr.  James  Bowers  [3]. 

Due  to  the  large  number  of  active  elements  in  IC  operational 
amplfiers,  a  device-level  computer  simulation  would  be  quite  expensive 
and  time  consuming.  Simulation  of  circuits  containing  several  op  amps 
would  not  be  a  reasonable  endeavor  using  a  model  of  this  order  of  com¬ 
plexity. 

A  macromodel  is  therefore  desired  which,  while  accurately  predicting 
the  performance  of  the  op  amp  under  all  circuit  conditions,  will  have 
a  minimum  level  of  complexity,  so  that  larger  circuits  containing  several 
op  amps  may  be  efficiently  simulated  in  reasonable  amounts  of  computer 
time. 

This  model  should  contain  all  of  these  common  op  amp  characteristics 

1.  Input  impedance,  bias  current,  offset  current  and  voltage. 

2.  Output  impedance  (DC  and  AC),  voltage  and  current  limiting. 

3.  Differential  and  common  mode  gain  and  rolloff  versus  frequency. 

4.  Non-symmetrical  slew  rate  and  step  response  feedthrough. 

The  model  should  be  applicable  to  both  standard  bipolar  input  op 
amps  as  well  as  devices  employing  field  effect  transistors  in  their  input 
stages,  and  should  accept  external  compensation  if  required. 

In  this  report,  a  model  is  developed  by  simplifying  a  previously 
suggested  slewing  model  to  a  two-diode  equivalent  circuit.  This  circuit 
is  then  modified  to  include  the  desired  input,  output,  and  differential 
features.  The  device  parameters  (frequency  response,  etc.)  is  derived 
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for  this  general  model,  then  the  model  parameters  are  solved  for  in 
terms  of  the  device  parameters.  Finally,  models  are  developed  for 
four  real  op  amps  and  extensively  tested. 

Volume  II  B:  "Macromodeling  of  Digital  Circuits,"  by  Ronald  Vogelsong 
and  Dr.  James  Bowers  [4]. 

A  NAND  gate  macro-model  is  developed  using  SPICE  which  allows  for 
the  complete  simulation  of  the  operation  of  the  logic  gate.  The  model 
derivation  requires  only  external  DC  voltage  and  current  measurements 
and  logic  delay  times,  but  the  model  developed  is  shown  to  fully  simulate 
the  high-frequency  response  of  the  device  as  well,  including  the  response 
of  the  device  to  interference  in  the  l-100MHz  range. 

In  addition  to  the  fully-qualified  model  developed  for  SPICE,  a 
simpler  model  is  constructed  for  use  in  the  SUPER*SCEPTRE  simulation  pro¬ 
gram  which  effectively  models  the  device  from  DC  up  to  normal  switching 
speed  limits. 

Finally,  a  look  is  taken  at  developing  a  SPICE  model  of  a  device 
constructed  from  a  large  number  of  logic  elements  fabricated  on  a 
single  chip",  a  TTL  5485  four  bit  comparator.  In  this  case,  the  overall 
complexity  is  such  that  a  modeling  procedure  is  developed  based  on  mini¬ 
mizing  the  complexity  and  hence  computer  (CPU)  time  of  the  circuit 
while  still  maintaining  the  logical,  input,  output,  and  propagation 
delay  characteristics. 


Volume  III  A:  Probabilistic  Analysis  of  Combinational  Circuits 

with  Random  Delays"  by  Abner  Ephrath  and  Dr.  Donald  Weiner  [5] 

Random  propagation  delays  are  encountered  in  digital  integrated  cir¬ 
cuits  due  to  fluctuations  in  the  fabrication  process.  These  delays  can 
be  further  increased  due  to  the  presence  of  electromagnetic  interference. 
System  performance  can  be  evaluated  from  the  expected  values  of  the  output 
signals.  Analytical  methods  for  determining  the  output  expected  values 
of  combinational  circuits  with  random  delays  are  developed  in  this  report. 

Given  the  input  expectations,  the  network  logic  functions,  and  p.d.f.'s 
of  the  delays  associated  with  the  gates  in  the  network,  it  is  shown  how  to 
obtain  the  output  expected  values.  Two  types  of  delay  elements  are  con¬ 
sidered:  1)  the  pure  delay  element,  whose  output  is  a  delayed,  but  un¬ 

distorted,  replica  of  the  input  and  2)  the  discriminating  delay  element, 
where  input  rise  and  fall  transitions  experience  different  delays.  Two 
degrees  of  network  complexity  are  dealt  with:  1)  tree-like  networks, 
in  which  there  is  only  one  path  from  every  network  input  to  any  network 
output  and  2)  networks  with  reconvergent  fanouts,  where  more  than  one 
path  exists  from  some  inputs  to  some  outputs. 

To  simplify  analysis  of  very  large  circuits,  an  approximate  model 
is  proposed  where  the  circuits  are  subdivided  into  large  logic  blocks. 

The  analytical  techniques  previously  derived  for  individual  gates  are 
then  applicable.  Various  strategies  for  characterizing  the  delays  of 
the  large  logic  blocks  are  considered  and  examined  by  means  of  computer 


simulations. 


Volume  III  B:  "Computer  Simulation  of  EMI  Effects  in  a  Differential 

Line  Receiver"  by  Zino  Chair,  Tushar  Dave,  Hamza  Ouibrahim, 
and  Dr.  Donald  D.  Weiner  [6], 

It  is  frequently  necessary  to  transmit  baseband  digital  data  over 
long  transmission  lines  such  as  system  interconnect  cables.  When  the 
cables  are  located  in  a  high  EMI  environment,  ordinary  integrated  logic 
circuits  cannot  be  used  to  synthesize  the  line  driver  and  receiver  because 
such  circuits  are  unable  to  discriminate  between  a  valid  signal  and  the 
externally  induced  interference.  One  solution  for  combating  the  inter¬ 
ference  is  to  convert  the  ground  referenced  digital  data  at  the  driver 
into  a  differential  signal  which  is  transmitted  over  a  pair  of  wires 
such  as  a  twisted-pair  line.  A  differential  line  receiver  is  used  at 
the  other  end.  Interference  which  appears  equally  on  each  wire  is 
referred  to  as  common-mode  interference.  Ideally,  the  line  receiver 
responds  only  to  the  differential  signal  transmitted  over  the  line  and 
rejects  the  common-mode  interference. 

This  report  describes  some  preliminary  results  dealing  with  the 
computer  simulation  of  EMI  effects  due  to  sinusoidal  interference  in  a 
DS  55115  differential  line  receiver.  Attention  is  focused  on  the  receiver 
as  opposed  to  the  line  driver,  because  the  receiver  is  significantly 
more  susceptible.  Consequently,  the  susceptibility  of  the  driver-receiver 
pair,  in  most  cases,  is  determined  by  that  of  the  receiver. 

The  computer  program  selected  for  the  investigation  was  SPICE 
(Simulation  Program  with  Integrated  Circuit  Emphasis)  which  was  exercised 
in  its  transient  mode.  Thus,  a  time-domain  analysis  was  performed  and 


time  waveforms  were  generated  at  all  nodes  of  interest.  Parameter 
values  of  the  semiconductor  devices  embedded  in  the  integrated  circuit 
for  the  receiver  are  not  published  in  the  interface  handbooks.  Reason¬ 
able  values  were  determined  by  adjusting  nominal  values  until  the  computer 
simulated  characteristics  closely  agreed  with  typical  performance  charac¬ 
teristics  supplied  by  the  manufacturer.  These  included  curves  of  1)  input 
current  vs.  input  voltage,  2)  output  voltage  vs.  common-mode  input  voltage 
3)  high  level  output  voltage  vs.  high  level  output  current,  4)  low  level 
output  voltage  vs.  low  level  output  current,  5)  output  voltage  vs.  dif¬ 
ferential  input  voltage,  6)  output  voltage  vs.  strobe  input  voltage, 
and  7)  power  supply  current  vs.  power  supply  voltage. 

The  sinusoidal  interference  was  impressed  at  the  two  data  inputs, 
the  strobe  terminal,  and  the  response-time  control  input.  Waveform  dis¬ 
tortion  in  the  output  signal  was  noted  for  suitably  strong  interferers. 
Waveform  parameters  for  describing  the  distortion  are  defined  and  plotted 
as  a  function  of  the  frequency  and/or  amplitude  of  the  sinusoidal  inter- 
ferer. 

Particular  attention  is  devoted  to  dc  offsets  caused  by  rectifica¬ 
tion  in  the  base-to-emitter  and  base-to-collector  transistor  junctions. 

The  input  voltage  at  which  the  output  changes  state  is  defined  as  the 
input  threshold  voltage.  The  rectification  of  RF  energy  causes  the  input 
threshold  voltage  to  change.  The  concept  of  using  this  change  to  define 
receiver  susceptibility  is  examined. 


In  addition  to  decreasing  the  noise  margin,  a  change  in  the 
input  threshold  voltage  can  cause  time  variations  in  the  receiver 
output  (i.e. ,  pulses  may  appear  shifted  in  time  and  some  pulses  may 
appear  longer  or  shorter).  The  resulting  jitter  constrains  the  minimum 
pulse  width  and,  therefore,  the  maximum  data  rate  which  can  be  used. 

This  effect  is  examined  as  a  function  of  the  amplitude  and/or  frequency 
of  the  sinusoidal  interference. 

Volume  IV  A:  "Coupling  of  Electromagnetic  Fields  into  Transmission 
Lines:  A  comparison  of  the  Transmission  Line  Model 
and  the  Method  of  Moments"  by  Robert  Abraham  and  Dr.  Clayton 
Paul  [7], 

The  coupling  of  electromagnetic  fields  onto  transmission  lines  is 
investigated.  The  transmission  line  model  with  distributed  sources 
is  employed  as  a  computationally  efficient  method  of  predicting  currents 
induced  on  transmission  wires  by  an  incident  field.  Results  obtained 
from  the  transmission  line  model  solution  are  compared  with  predictions 
made  by  the  more  rigorous,  but  much  less  efficient,  method  of  moments 
technique.  Two  user-oriented  computer  codes,  using  different  expansion 
and  testing  functions,  were  selected  to  provide  the  method  of  moments 
solution.  Both  the  prediction  accuracy  and  limitations  of  the  trans¬ 
mission  line  model  are  explored  in  depth  using  a  carefully  selected 
transmission  line  structure.  The  structure  is  modified  slightly  to 
illustrate  several  important  characteristics  of  the  transmission  line 


model  solution.  Low  frequency  limitations  of  the  method  of  moments 
solution  are  investigated.  In  addition,  the  differences  between  the 
two  method  of  moments  formulations  are  found  to  have  a  significant 
effect  on  the  integrity  of  their  individual  solutions.  These  differences 
are  illustrated.  Finally,  the  practical  use  of  the  transmission  line 
model  as  an  effective  method  for  predicting  the  coupling  of  an  elec¬ 
tromagnetic  field  into  terminal  devices  is  discussed. 

Volume  IV  B:  "Shield,  A  Digital  Computer  Program  for  Computing  Cross¬ 
talk  Between  Shielded  Cables"  by  Dr.  Clayton  Paul  [8]. 

This  report  contains  the  description  and  verification  of  a  digital 
computer  program,  SHIELD,  to  be  used  in  the  prediction  of  crosstalk 
in  transmission  lines  consisting  of  unshielded  wires  and/or  shielded 
cables.  The  line  may  be  above  a  ground  plane  (Type  1)  or  within  an 
overall,  circular,  cylindrical  shield  which  may  be  solid  or  braided  and 
a  wire  (the  shielded  wire)  located  concentrically  on  the  axis  of  the 
shield.  All  wires  may  be  stranded  and  all  conductors  are  treated  as 
imperfect  conductors;  that  is,  their  per-unit-length  impedances  are  non¬ 
zero.  Through-braid  coupling  for  braided  shields  as  well  as  diffusion 
for  both  types  are  included  in  the  model.  The  shielded  cables  may  have 
exposed  sections  at  either  end  (pigtail  sections)  in  which  the  shielded 
wire  is  not  covered  by  the  shield.  Over  these  pigtail  sections,  a  pig¬ 
tail  wire,  parallel  to  the  shielded  wire,  connects  the  shield  to  the 
reference  conductor  at  that  end  via  either  a  short  circuit  or  an  open 
circuit.  These  pigtail  sections  are  included  in  the  representation  to 


simulate  the  common  practice  of  terminating  a  shielded  cable  in  a 
connector  via  these  pigtail  wires.  The  pigtail  sections  may  be  of 
different  lengths.  The  program  is  written  in  FORTRAN  IV  and  should 
be  implementable  on  a  wide  range  of  digital  computers. 

Volume  IV  C:  "Crosstalk  in  Twisted  Wire  Circuits"  by 
Marty  Jolly  and  Dr.  Clayton  Paul  [9], 

An  investigation  of  wire-to-wire  crosstalk  is  presented  in  an  attempt 
to  accurately  predict  the  magnitude  of  the  voltage  coupled  to  a  twisted 
wire  pair.  A  sensitivity  analysis  of  the  twisted  wire  pair  is  conducted 
to  determine  the  effect  that  line  twist  has  on  the  coupled  differential 
mode  voltage  when  the  twisted  pair  is  connected  in  an  unbalanced  configur¬ 
ation.  An  improved  computer  model  is  developed  to  aid  in  the  prediction 
of  the  voltage  coupled  to  the  twisted  wire  pair  when  it  is  terminated 
in  low  impedance  loads.  Another  circuit  configuration  which  is  investig¬ 
ated  consists  of  a  differential  line  driver  connected  to  a  differential 
line  receiver  by  a  twisted  wire  pair.  Predictions  of  the  common  mode 
voltage  coupled  to  the  twisted  pair  are  computed  with  a  single  wire  model 
and  a  comparison  is  made  with  experimental  results. 

Volume  IV  D:  "Modeling  Crosstalk  in  Balanced  Twisted  Pairs"  by 
Dawn  Koopman  and  Dr.  Clayton  Paul  [10] 

The  coupling  of  electromagnetic  fields  to  twisted  wire  pair  trans¬ 
mission  lines  is  investigated.  A  simple  low-frequency  model  to  predict 
the  crosstalk  of  a  three-conductor  line  is  described.  This  low-frequency 
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model  is  used  to  demonstrate  how  the  twisted  pair  reduces  crosstalk. 

The  sensitivity  of  the  crosstalk  from  the  single  wire  to  unbalanced  twisted 
pair  is  examined.  The  crosstalk  from  the  unbalanced  twisted  pair  to 
unbalanced  twisted  pair  is  measured  and  explained.  Finally,  the  cross¬ 
talk  in  the  single  wire  to  balanced  twisted  pair  is  measured  to  illustrate 
that  balancing  the  terminal  configurations  of  the  twisted  pair  can  further 
reduce  crosstalk. 

Volume  IV  E:  "Modeling  of  Transmission  Lines:  A  Comparison  of  Lumped- 
Circuit  Iterative  Models  and  the  Transmission  Line  Model" 
by  Woodrow  Everett  III  and  Dr.  Clayton  Paul  [11]. 

A  comparison  between  the  predictions  of  various  lumped-circuit  iter¬ 
ative  models  of  transmission  lines  and  that  of  the  transmission  line  model 
is  given  for  both  the  two  and  three-conductor  cases.  For  frequencies  such 
that  the  line  length  is  approximately  1/10  of  a  wavelength,  or  less  in 
one  section,  the  lumped-circuit  models  yield  predictions  which  are  within 
±3dB  of  the  transmission  line  model  predictions,  and  which  are  essentially 
independent  of  the  value  of  load  impedance.  For  frequencies  where  the  line 
length  is  greater  than  approximately  1/10  of  a  wavelength,  more  than  one 
section  of  the  lumped  models  must  be  used  to  yield  accurate  predictions, 
but  the  prediction  accuracy  may  be  significantly  affected  by  the  load 
impedance  value.  The  range  of  accurate  predictions  (in  terms  of  electri¬ 
cal  length)  may  be  increased  by  adding  more  sections  of  a  lumped-circuit 
model,  but  this  range  increases  in  a  nonlinear  fashion.  Because  of  this 
nonlinear  relationship,  it  is  difficult  to  find  a  satisfactory  method  for 


estimating  the  number  of  sections  of  a  particular  lumped-circuit  model 
that  must  be  used  in  order  to  obtain  a  desired  range  of  accurate  pre¬ 
dictions  . 

Volume  IV F:  "Prediction  of  Crosstalk  in  Flatpack  Coaxial  Cables"  by 
Wayne  Beech  and  Dr.  Clayton  Paul  [ 1 1 A] . 

The  electromagnetic  crosstalk  present  in  a  class  of  cables  known  as 
flatpack,  coaxial  cables  is  investigated.  The  multiconductor  trans¬ 
mission  line  equations  are  used  to  derive  a  model  of  a  general  flatpack, 
coaxial  cable  consisting  of  n  cables.  This  model  is  then  implemented  as  a 
digital  computer  program  to  allow  simulation  of  the  crosstalk  levels  present 
in  the  cable.  This  simulation  is  compared  to  experimental  results  to  prove 
that  the  model  is  valid  and  can  accurately  predict  the  crosstalk  levels. 

Also  the  effect  of  common  impedance  coupling,  the  presence  of  drain  wires 
to  allow  for  connection  to  the  shields  of  the  cable,  and  the  presence  of 
pigtail  sections  (exposed  sections  of  wire)  is  examined  to  determine  their 
effect  on  the  overall  electromagnetic  crosstalk  present  in  the  cable. 

Volume  V;  "Model  Revision  for  the  IEMCAP"  by  Dr.  Clayton  Paul  [12]. 

This  report  contains  recommendations  for  changes  in  the  wire-to-wire 
coupling  subroutine  (WTWTFR)  in  the  Intrasystem  Electromagnetic  Compati¬ 
bility  Analysis  Program  (IEMCAP).  The  changes  are  intended  to  accomplish 
four  objectives:  (1)  to  modularize  that  subroutine,  (2)  to  speed  up  com¬ 
putation,  (3)  to  provide  models  which  have  more  sound  theoretical  basis 
than  those  presently  included  in  IEMCAP,  and  (4)  to  correct  certain  errors 


presently  in  the  current  models  in  IEMCAP.  The  first  objective  -  to  modular¬ 
ize  the  subroutine  -  is  important  from  the  standpoint  of  future  maintenance 
of  the  code.  The  other  three  objectives  concern  relatively  accurate  and 
speedy  predictions.  It  is  important  that  the  models  have  a  sound  theoretical 
basis  rather  than  have  the  ability  to  predict  only  certain  limited,  empiri¬ 
cal  data.  One  then  has  some  confidence  that  the  models  will  predict  some, 
as  yet,  uninvestigated  situation. 

The  subroutine  which  was  initially  delivered  by  the  contractor  con¬ 
tained  numerous  theoretical  inconsistencies  and  model  prediction  errors. 

These  were  thought  to  have  been  corrected  via  a  completely  rewritten  sub¬ 
routine.  The  models  in  that  revised  subroutine  were  based  on  the  modeling 
efforts  to  that  date.  Shielded  wires  and  twisted  pairs  had  not  been 
extensively  investigated  from  the  standpoint  of  models  for  predicting 
crosstalk.  Thus  models  for  those  portions  of  the  subroutine  were  based  on 

limited  modeling  data  which  existed  at  the  time. 

Since  that  initial  revision,  several  modifications  were  made  to 
correct  coding  errors  and  to  add  features.  The  original 
revision  was  designed  to  handle  pigtails  on  shielded  wires  which  were 
three  inches  (3")  in  length.  No  provision  was  made  for  peripherally 
bonded  shields  (zero-length  pigtails)  or  variable  lengths  of  pigtails 
from  shield  to  shield.  Recently  a  revision  was  made  to  WTV1TFR  to  allow 
for  variable  length  pigtails.  It  has  been  determined  that  this  revision 
did  not  correctly  handle  pigtails  for  inductive  coupling.  In  the  original 
revision,  the  restriction  of  3"  pigtails  on  all  shield  terminations 
allowed  an  optimization  of  the  code  which  would  not  have  been  possible 
with  variable  length  pigtails.  Since  that  code  structure  and  the  models 
relied  heavily  on  the  restriction  that  all  shield  pigtail  terminations 
be  the  same,  it  is  not  a  simple  matter  to  change  that  code  to  now  handle 


variable  length  pigtails.  Moreover,  adding  the  capability  of  variable 
length  pigtails  by  modifying  the  original  code  would  not  take  advantage 
of  the  optimization  which  the  assumption  of  3"  pigtails  allowed.  It 
would  seem,  therefore,  that  one  should  rewrite  and  optimize  the  structure 
of  the  code  for  the  variable  length  pigtail  case  rather  than  patching 
up  the  original  revision. 

Much  additional  work  has  been  done  on  the  modeling  of  crosstalk 
involving  shielded  wires  and  twisted  pairs  since  that  original  code  re¬ 
vision.  This  additional  work  has  shown  that  some  of  the  models  for  these 
cases  contained  in  the  original  revision  are  not  correct  and  some  models 

were  unnecessarily  complicated  for  the  prediction  accuracies  which  one 
would  reasonably  expect  on  practical  systems.  In  addition,  some 
long-held,  fundamental  notions  concerning  the  superposition  of  induc¬ 
tive  and  capacitive  coupling  (which  were  fundamental  to  the  original 
revision)  were  shown  to  be  incorrect. 

Thus  because  of  the  additional  modeling  experience  and  the  discovery 
of  the  error  of  certain  fundamental  model  premises,  the  WTWTFR  models 
need  modification.  Because  of  the  extensiveness  of  these  modifications 
as  well  as  the  other  reasons  alluded  to  above,  it  seems  advisable  to 
completely  rewrite  the  WTWTFR  subroutine.  In  doing  so  it  is  also  advis¬ 
able  that  it  be  modularized  and  optimized.  Modularization  is  critically 
important  from  the  standpoint  of  future  maintenance.  It  is  possible, 
but  very  difficult,  for  anyone  who  is  not  intimately  (and  currently) 
familiar  with  the  present  code  to  make  changes  to  it  much  less  track 
the  effects  of  those  changes  through  the  code.  Modularization  would 
remedy  this  problem  to  a  large  degree. 
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Volume  VI:  "Basic  EMC  Technology  Advancements  for  CJ  Systems"  by 
Drs.  Clayton  Paul,  Donald  Weiner,  and  James  Whalen. 

This  research  effort  was  devoted  to  eliminating  some  of  the  basic 
technological  gaps  in  the  two  important  areas  of  1)  electromagentic 
effects  (EM)  on  microelectronic  circuits  and  2)  EM  coupling  and  testing 
The  results  are  presented  in  fourteen  reports  which  have  been  organized 
into  six  volumes.  The  reports  are  briefly  summarized  in  this  volume. 

In  addition,  an  experiment  is  described  which  was  performed 
to  demonstrate  the  feasibility  of  applying  several  of  the  results 
to  a  problem  involving  electromagnetic  interference.  Specifically, 
experimental  results  are  provided  for  the  randomness  associated  with 
1)  crosstalk  in  cable  harnesses  and  2)  demodulation  of  amplitude 
modulated  (AM) signals  in  operational  amplifiers.  These  results  are 
combined  to  predict  candidate  probability  density  functions  (pdf’s) 
for  the  amplitude  of  an  AM  interfering  signal  required  to  turn  on  a 
light  emitting  diode.  The  candidate  pdf’s  are  shown  to  be  statisti¬ 
cally  consistent  with  measured  data. 


II.  Theoretical  Discussion  of  Feasibility  Experiment 

Since  research  was  conducted  on  a  number  of  different  topics,  it  was 
decided  to  perform  an  experiment  that  demonstrated  the  feasibility  of  applying 
several  of  the  results  to  a  problem  involving  electromagnetic  interference 
(EMI).  The  research  areas  encompassed  by  the  experiment  were: 

1)  EMC  modeling  and  analysis  using  a  probabilistic  approach, 

2)  crosstalk  in  cable  harnesses, 

and 

3)  weakly  nonlinear  behavior  of  operational  amplifiers. 

This  chapter  discusses  theoretical  considerations  associated  with  the  experi¬ 
ment  . 

2 . 1  Probabilistic  Model  for  a  Typical  EMI  Problem 

A  typical  EMI  problem  is  illustrated  in  Fig.  2.1-1. 


Fig.  2.1-1.  A  typical  EMI  problem 

The  desired  signal  for  the  equipment  is  denoted  by  d(t).  It  is  transmitted  to 
the  equipment  via  coupling  path,  P^,  and  is  received  as  the  signal,  s(t).  On 


the  other  hand,  the  interfering  signal,  j(t),  is  coupled  to  the  equipment  via 
path,  P?  and  arrives  as  the  signal,  i(t).  Uncertainties  and/or  randomness  in 


d(t),  j(t),  and  result  in  a  random  behavior  for  both  s(t)  and  i(t). 

Also,  variability  in  such  things  as  electrical  components  and  the  positioning 
of  wires  and  cables  result  in  a  variation  of  performance  from  one  piece  of 
"identical"  equipment  to  the  next. 

In  the  probabilistic  approach  [14,15,161,  s(t)  and  i(t)  are  both  treated 
as  sample  functions  from  random  processes.  In  addition,  the  equipment  is 
viewed  as  a  sample  from  an  ensemble  of  randomly  distributed  "identical"  equip¬ 
ments  whose  EMI  performance  is  characterized  in  terms  of  a  random  susceptibil¬ 
ity  level.  The  problem  is  then  analyzed  by  determining  the  probability  of 
EMI  [16]. 

2 . 2  Description  of  the  Experiment 

A  block  diagram  of  the  experiment  is  shown  in  Fig.  2.2-1.  With  respect 
to  the  discussion  of  section  2.1,  the  equipment  interfered  with  is  the  printed 
circuit  board  (PCB) .  The  network  in  the  PCB  was  designed  such  that  the  light 
emitting  diode  (LED)  turned  on  when  the  input  voltage,  V^,  exceeded  the  speci¬ 
fied  threshold,  The  three  stages  of  amplification  served  to  amplify  the 

input  signal  so  that  the  output  voltage,  V^»  was  just  large  enough  to  light 
the  LED  when 

Normally,  the  LED  was  intended  to  be  off.  During  the  experiment  an  inter 
fering  signal  was  coupled  to  the  PCB  through  a  25-wire  random  cable  bundle. 

EMI  occurred  when,  for  large  enough  amplitudes  of  the  interfering  signal,  the 
LED  was  turned  on.  Since  the  LED  was  normally  off  during  the  experiment,  the 
desired  signal  in  Fig.  2.1-1  was  identically  zero.  Specifically, 

d(t)  =  s(t)  =  0.  (2.2-1) 


The  interfering  signal  was  selected  to  be  the  amplitude  modulated  (AM)  signal 
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with  sinusoidal  modulation  given  by 


j(t)  =  A  [1  +  m  cos  (2frf  t)]  cos  27Tf  t.  (2.2-2) 

i  me 

The  numerical  values  of  the  parameters  in  Eq.  (2.2-2)  used  in  the  experiment 

were  m  =  0.5,  f  =  1kHz,  and  f  =  4MHz  and  15  MHz. 
m  c 

The  circuit  schematic  of  the  PCB  is  shown  in  Fig.  2.2-2.  All  three 

stages  employ  operational  amplifiers.  The  first  two  stages  are  inverting 

amplifiers  with  DC  gains  of  approximately  10  and  100,  respectively.  The  third 

stage  is  a  noninverting  buffer  with  a  DC  gain  of  approximately  10.  As  a 

result,  the  overall  DC  voltage  gain  from  the  input  voltage  V to  the  output 

voltage  is  approximately  10,000. 

The  AM  interfering  signal  contains  frequency  components  at  f  -  f  , 

c  m 

f  ,  and  f  +  f  .  Since  f  <<  f  ,  all  of  these  frequencies  are  either  at  or 
c  c  m  me 

very  close  to  the  carrier  frequency.  The  low  pass  filter  following  the  first 
stage,  which  consists  of  and  C^,  has  a  3-dB  cut-off  frequency  of  approxi¬ 
mately  1700  Hz.  Therefore,  if  the  circuit  behaved  in  a  linear  fashion,  the 
interfering  signal  would  be  severely  attenuated  at  the  output  of  the  PCB.  As 
a  result,  the  interference  would  not  light  the  LED. 

However,  operational  amplifiers  are  weakly  nonlinear  in  the  amplification 
region.  Therefore,  during  the  experiment,  the  first  stage  served  to  demodu¬ 
late  the  interfering  AM  signal.  The  DC  component  of  the  modulation  was  pre¬ 
vented  from  appearing  at  the  output  by  the  blocking  capacitor,  C^.  On  the 
other  hand,  the  circuit  amplified  the  1kHz  sinusoidal  component  of  the  modu¬ 
lation  which  fell  well  within  the  passband  of  the  three  stages.  The  LED  was 
observed  to  light  for  large  enough  amplitudes  of  the  interfering  signal. 
Lighting  of  the  LED  depended  on  two  distinct  phenomena: 


1)  the  interfering  signal  was  coupled  to  the  PCB  by  means  of  crosstalk  which 
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occurred  within  the  cable  bundle  and  2)  a  1  kHz  sinusoidal  component  was  gen¬ 


erated  due  to  the  weakly  nonlinear  behavior  of  the  operational  amplifier  in 
the  first  stage.  In  chapter  III  it  is  shown  that  changing  the  wrapping  of 
the  cable  bundle  changes  the  strength  of  the  interfering  signal  coupled  to 
the  PCB.  In  Chapter  IV  it  is  shown  that  changing  the  operational  amplifier 
changes  the  amplitude  of  the  generated  sinusoidal  component  at  1  kHz.  The 
two  changes  together  result  in  various  fluctuations  in  the  amplitude  of  the 
interfering  signal  at  needed  to  light  the  LED. 

Separate  experiments  were  carried  out  to  determine  the  probability  den¬ 
sity  functions  of  the  quantities  that  characterize  1)  crosstalk  within  the 
cable  bundle  and  2)  demodulation  of  the  interfering  signal.  These  density 
functions  were  then  combined  analytically  to  determine  the  probability  den¬ 
sity  function  of  the  voltage  at  which  results  in  a  specified  threshold  vol¬ 
tage  at  Vy 

Another  separate  experiment  was  conducted  in  which  the  cable  wrappings 
and  operational  amplifier  in  the  first  stage  were  changed  many  times.  Cable 
wrappings  and  operational  amplifiers  from  the  previous  experiments  were  not 
used.  The  voltages  at  required  to  produce  a  specified  threshold  voltage  at 
were  recorded.  Finally,  a  statistical  test  was  employed  to  determine 
whether  the  data  for  were  statistically  consistent  with  the  analytically 
obtained  probability  density  function  for  . 


2.3  Analysis  of  Weakly  Nonlinear  Effects  in  the  Operational  Amplifier 

The  interfering  signal  at  the  input  to  the  cable  bundle  is  an  AM  waveform 
of  the  form 


j(t)  =  A  [  1+m  cos  ( 2tt f  t)]  cos  2tt f  t 


(2.3-1) 


Expansion  of  Eq. (2.3-1)  and  application  of  a  well-known  trigonometric 
identity  results  in 

mA1 

j(t)  =  A  cos  (2iTf  t)  +  cos  ( 2tt ( f  -f  )t]  (2.3-2) 

i  c  z  cm 

mA 

+  — 7T-  cos  [  27T  ( f  +f  )t]. 
z  cm 


Let  the  cross  talk  within  the  cable  bundle,  from  to  V^,  be  characterized 


by  the  linear  voltage  transfer  function,  G(f).  Since  f  =  1kHz  while 

m 


f  =  4mHz  or  15  mHz, 
c 


f  -  f  s  f  +  f  s  f  . 
c  m  c  m  c 


(2.3-3) 


Therefore,  it  is  reasonable  to  assume 


j° 


G(f c_fm>  a  G<fc  +  fm>  ?  G(fc>  "  lG(fc)le  * 


(2.3-4) 


Let  the  amplitude  of  the  sinusoidal  carrier  at  the  output  of  the  cable  bundle 


be  denoted  by  k^.  It  follows  that 


A2  =  |G(fc) |  Aj .  (2.3-5) 

Because  the  three  components  of  Eq. (2.3-2)  experience  approximately  the 
same  gain  and  phase  shift,  it  follows  that  the  output  of  the  cable  bundle  is 
a  sinusoidally  modulated  AM  signal  given  by 


i(t)  =  A„  [  1  +  m  cos(2Tif  t)]  cos(2Tif  t  +  o  ) 
t-  m  c  c 


mA„ 


A„  cos(2Ttf  t  +  o  )  +  — —  cos[2ir(f  -f  )t  +0  ]  (2.3-6) 
2  c  c  2  cm  c 


mA„ 


+  — —  cos[2TT(f  +f  )t  +o  ] 


To  simplify  Eq. (2.3-6),  it  is  convenient  to  define 


f  =  f  -  f  ,  E  = 
1  c  m  l 


jo  mA  jo 

l*,|e  ' 


f2  '  fc  ’  E2 


f  =  f  +  f  ,  E.  = 
3  c  m  3 


JO  3° 

|e2|b  c  ■  A2  e  c 


J°c  »A2  j°c 
|e3|b  -  —  e 


(2.3-7) 


Eq.  (2.3-6)  may  then  be  expressed  as 


i(t)  =  Y  |E  |cos(2tt  f  t  +  O  ) 
q=l  q  q 


,  3  j2TTfqt 

'  2  2  ,  Eq  6 

q=-3 
q^O 


(2.3-8) 


where 


-J°, 


E  =  E  e 
-q  1  q1 


and  f  =  -  f  . 
-q  q 


(2.3-9) 


The  operational  amplifiers  in  the  second  and  third  stages  of  the  PCB 
are  assumed  to  behave  linearly.  However,  the  weakly  nonlinear  operational 
amplifier  in  the  first  stage  is  assumed  to  possess  a  second-order  nonlinear¬ 


ity  (i.e.,  the  input  signal  driving  the  first  stage  is  assumed  to  be  small 
enough  such  that  higher-order  nonlinearities  may  be  considered  to  be 
negligible.)  In  the  experiment  it  is  the  second-order  component  at  fre¬ 
quency,  f  ,  which  is  of  interest. 

From  Eq. (2.3-8)  it  is  obvious  that  i(t)  contains  the  six  input  frequen¬ 


cies  f_3  =  -f3>  f_2  =  -f2,  f_L  =  -fj,  fj,  f2,  f3.  Let 


denote  the  number 


th 


of  times  the  frequency  f  appears  in  a  particular  n  -order  frequency  mix. 


m 


All  possible  n  -order  frequency  mixes  are  then  represented  by  the  frequency 


mix  vector 


in  —  (m  3>  m  »  1^2  >  m^) 


(2.3-10) 


where 


1  \  =  n- 
k=-3  “ 


(2.3-11) 


The  intermodulation  frequency  corresponding  to  m  is 


3 

£  -  I 


\fk  =  ^ml  “  m-Pfl  +  (m2_m-2^f2  +  (m3-m_3)f3.  (2.3-12) 


In  terms  of  the  nonlinear  transfer  function  approach  [17],  the  n  -  order 
response  associated  with  m  to  the  input  given  by  Eq.  (2.3-8)  is 

(n;m)  (m  +m  )  (m  +m  ~)  (m.+m  _) 

*»<'•=>  -^1  lEl'  ^  3 


I H  (m)  I  cos  [2irf  t  +  O  +  ip  (m)  ] 
n  —  '  m  m  n  — 


(2.3-13) 


where  (n;m)  is  the  multinomial  coefficient  defined  by 


(n;m)  = 


(m_3 ! ) (m_2 1 ) (m_x ! ) (m1 ! ) (m2 ! ) (m3 ! ) 


Hn(m)  is  the  n  -order  nonlinear  transfer  function  given  by 


V-)  =  lHn^>|  e 


j'^n(m) 


(2.3-14) 


(2.3-15) 


=  Hn(f_3,...  ,f_3,...,  f_i***»f  _!»•••» f  !»•••»  fj» *  f  •  *  .  f 


.V'-’  -*•  «V >•>  -V  W/?/  VA- VV vWV 


°m  =  ki_3  Vc  =  I(mrn,-1>+  (ra2-m_2)+  (m3"m-3)]0c- 


(2.3-16) 


Eqs.  (2.3-13)  through  (2.3-16)  are  now  specialized  to  the  second-order 
response  at  frequency,  f^. 

With  respect  to  the  frequencies  defined  in  Eq.  (2.3-7),  the  frequency 
mixes  which  result  in  a  second-order  reponse  at  fm  are: 

a) (f2-fL) 

b) (f3-f2). 


The  frequency  mix  vector  corresponding  to  the  mix  (f2-f^) 

m^  =  (0,0, l, 0,1,0) 

while  that  corresponding  to  the  mix  (f^  -  f^)  is 

^2  =  (0,1, 0,0, 0,1). 

The  associated  multinomial  coefficients  are 


(2;mL)  =  (2;m2)  =  2. 


(2.3-17) 


(2.3-18) 


(2.3-19) 


In  addition. 


0  =  (m_  -  m  , )  o  =0 

m^  Z  —1  c 


a  =  (m  -  m  „)o  =  0. 

E-2  3  “2  c 


It  follows  that  the  total  second-order  response  at  f  is 

m 


(2.3-20) 


y2(t;fm)  =  y2(t;H1)  +  y2(t;H2) 


(2.3-21) 


=  |E1 1  |E2  |  |  H2  (-f  1  ,f2)  j  cos  [  2tt  (f  2-f  x )  t  +  \p2  (m^] 


+  |E2  |  |E  j  | H2 (— f 2 , f  )  j  cos [ 2tt ( f 3— f 2 ) t  +  t^2  ( rn 2 )  ] 


Finally,  use  of  Eqs.  (2.3-7)  in  Eq. (2.3-21)  yields 

2 


mA„ 


y„(t;f  )  =  —7T 
2  m  2 


mA,. 


1 H2 (— f  i » f  2)  1  C°S  ^2lTfmt  +^2^— 1^ 


~2~~  iH2(-_f2’f3^  COS  ^27Tfmt  +^2(‘-2^’ 


(2.3-22) 


The  above  result  can  be  simplified  even  further. 
Because  of  Eq.  (2.3-3), 


f  =  f_  s  f_  =  f  . 
1  3  2  c 


(2.3-23) 


Consequently , 


H2(~fl,f2)  2  H2(_f2*f3)  S  H2(-fc‘fc)* 


(2.3-24) 


(-f  ,  f^)  is  associated  with  the  frequency  mix,  (f c~f  ) ,  which  yields 
a  DC  component.  Therefore,  ^(-f^f  )  3  pure^y  real  quantity.  Since 

H2 (— f  i » f 2 ^  and  H2^~^2’^3^  are  *30t^  approximately  real, 


(2.3-25) 


Therefore,  Eq. (2.3-22)  simplifies  to 


y0(t;f  )  ~  mA„  !  H„  (— f  ,f  )|  cos  2nf  t. 
2  m  2  2  ^  0  m 


c  c 


m 


(2.3-26) 


Utilizing  Eq.  (2.3-5),  the  amplitude  of  the  sinusoidal  component  at  f  is 

m 


given  by 

,2 


"^2  iV-VVi  =miG^fc)|2  Al2lH2("fc’fc) 


(2.3-27) 


2 . 4  Probabilistic  Analysis  of  Experiment 

Instead  of  analyzing  the  circuit  in  Fig.  2.2-2  stage  by  stage,  it  is 
convenient  to  treat  the  entire  PCB  as  a  single  system  having  input  and 

29 


output  V^.  This  is  readily  accomplished  by  defining  the  second-order  nonlinear 
transfer  function  of  the  previous  section  such  that  it  relates  the  AM  inter¬ 
ference  at  V„  to  the  1kHz  tone  at  V„.  In  this  way,  H„(-f  ,f  )  accounts  for 
z  J  zee 

the  linear  amplification  provided  by  the  second  and  third  stages  in  addition 
to  the  weakly  nonlinear  behavior  of  the  first  stage.  With  this  definition, 
y^Ctjf^)  of  the  previous  section  corresponds  to  V^(t).  I  :t  denote  the 
amplitude  of  the  1kHz  tone  at  V^.  Then  Eq.  (2.3-26)  can  be  written  as 

V_(t)  =  A_  cos  2-rrf  t  (2.4-1) 

dim 

where 

A~  =  mA  2 [ H  (-f  ,f  ) | 

J  z  1  2  c  c  1 

=  m|G(fc)|2  A12|H2(-fc,fc)|.  (2*A_2) 

During  the  experiment  was  increased  until  the  LED  was  lit.  Denote  the 
values  of  A^  and  A^  at  which  the  LED  switched  on  by  A^  and  A^,  respectively. 
From  Eq.  (2.4-2),  and  A^  are  related  by 

A3T  =  mlGM2  AlT2|H2("fc»fc)|*  (2'4_3) 


The  experiment  was  repeated  many  times.  For  each  run,  the  cable  was 
rewrapped  and  a  different  operational  amplifier  was  placed  in  the  first  stage. 
The  rewrappings  resulted  in  random  values  for  | G (f  ) | .  Similarly,  changing 
the  operational  amplifier  yielded  random  values  for  | H2 (-f  , fc) | .  In  this 
section  the  probability  density  function  for  A^  is  derived  given  a  specified 


value  for  A_  and  the  probability  density  functions  for  | G ( f  )|  and 

Jl  ^ 


M-f  ) 


From  Eq.  (2.4-3),  is  given  by 


IT  m 


|G(fc)|[|H2(-fc,fc)|]2 


(2.4-4) 


For  ease  of  analysis,  let  the  random  variables  associated  with  ]  G ( f c)  1  an(* 

| (— f c » f c ) 1  '3e  denoted  by  G  and  H,  respectively.  Also,  let  values  assumed 
by  the  random  variables  A^,  G,  and  H  be  denoted  by  a^,  g,  and  h  ,  respective 
ly.  The  probability  density  function  of  A^  is  obtained  by  considering  the 
following  two  functions  of  two  variables 


Alx  =  k(G,H)  = 


W  =  £ (G,H)  =  H. 


(2.4.5) 


In  terms  of  the  values  assumed  by  the  random  variables, 

A3T  'i  | 

aIT  '  k(s'h)  "  IT  — tV 


w  =  £(g,  tj>  =  k 


The  Jacobian  of  the  transformation  in  Eq.  (2.4-6)  is 


(2.4-6) 


J(g»h)  = 


:3k(g,h)  )k(g,h) 

?g  r.k 

’«  (g>h)  -  ■  (g,h) 

:<g  3h 


3T  1 


g'  [h] 


•  rA 

I  _  i 

,  2  m 


A  ^ 
3T  ^ 

m 


(2.4-7) 


Let  the  probability  density  functions  for  G  and  H  be  denoted  by  f  (g) 

Cj 


and  f  (h) ,  respectively.  Also,  let  the  joint  probability  density  function 
H 


for  the  random  variables  A.„  and  W  be  denoted  by  f.  (a  ,w)  .  Solution  of 

IT  A1tW  IT 


Eqs.  (2.4-6)  for  g  and  h  in  terms  of  a^  and  w  yields 


g  = 
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a^rjJw] 


h  =  w. 


(2.4-8) 


Assuming  G  and  H  to  be  statistically  independent,  it  can  be  shown  that  [18] 


£A1TU  (aLC“>  ■ 


|  fG(^fH(h> 

A3T 

g  = 

-  01 

1 


(2.4-9) 
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h  =  W 
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3T 


alT[w] 


Finally,  f  (a  )  is  obtained  by  integrating  f  (a  ,w)  over  all  possible 

1  i-  I L 


values  of  w.  This  results  in 

l. 


fA  (a  )  = 

V  1T 


3T 
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IT 


fw]  2 


X  fr 


'h 


3T 


a1T[w] 


r  fH(w)dw.  (2.4-10) 


h\  h 


liccause  w  =  h,  a  more  convenient  form  of  Eq.  (2.4-10)  is 


AlT(aLT) 


A 


3T 


L_  f  _L_  f 

J-UO  ^  ^ 


alT2 


[h] 


3T 


m 


■»  i  ) 

- — r  jfu(h)dh .  (2.4-11) 


alT[h] 


It  is  seen  that  knowledge  of  A3T  and  the  probability  density  functions 


for  G  and  H  is  sufficient  to  determine  the  probability  density  function  for 


Crosstalk  Variations  in  Cable  Harnesses 


Inadvertent  electromagnetic  coupling  between  wires  in  cable  bundles 
(crosstalk)  is  recognized  as  a  potentially  serious  cause  of  degradation  in  the 
performance  of  an  electronic  system.  Numerous  methods  of  combatting  cross¬ 
talk  have  been  utilized  in  the  past  [191.  Those  methods  involve  the  use  of 
shielded  and/or  twisted  wires  as  well  as  filters.  In  many  cases  a  simple 
rerouting  of  the  wires  into  other  cable  harnesses  can  bring  about  crosstalk 
reduction.  Of  course,  the  degree  to  which  rerouting  wires  reduces  the  cross¬ 
talk  depends  on  the  spectral  content  of  the  signals  on  adjacent  wires  in  the 
harness  as  well  as  the  degree  of  coupling  between  those  wires. 

In  recent  years,  methods  of  mathematically  modeling  this  crosstalk  to 
determine  its  degree  of  severity  have  been  developed  and  studied.  All  of 
these  modeling  techniques  rely  on  the  theory  of  the  distributed  parameter, 
multiconductor  transmission  line  model.  For  frequencies  such  that  the  cable 
length  is  sufficiently  short,  electrically,  lumped  parameter  approximations  to 
this  basic  model  are  sufficient.  For  controlled  configurations  wherein  the 
relative  wire  spacings  are  known  and  controlled  along  the  cable  bundle,  this 
model  has  been  shown  to  yield  accurate  predictions  of  the  crosstalk  [  19]  . 

Although  current  harness  technology  is  making  more  frequent  use  of  the 
controlled  characteristic  flatpack  types,  there  exist  a  large  proportion  of 
the  cable  harnesses  which  may  be  classified  as  random  cable  bundles.  In 
these  traditional  types  of  cable  harnesses,  the  wires  are  held  in  close  proxi¬ 
mity  by  lacing  cord,  cable  ties,  etc.  The  exact,  relative  positions  of  the 
wires  are  unknown  and  moreover  vary  along  the  cable  bundle.  Two  supposedly 
identical  harnesses  in  two  supposedly  identical  systems  will,  in  reality  ex¬ 


hibit  quite  different  crosstalk  properties 


Also,  these  crosstalk  properties  depend  quite  strongly  on  the  terminal 
loads  attached  to  the  ends  of  the  wires.  In  a  recent  study  [20],  this  effect 
was  shown  to  be  significant.  A  25-wire  random  cable  bundle  shown  in  Fig.  3-1 
was  constructed  of  #22  gauge  stranded  wires  suspended  2.4  cm  above  a  ground 
plane.  The  cable  length  was  4  m  and  the  wires  were  inserted  into  connectors 
at  each  end  as  shown  in  Fig.  3-2.  These  connectors  were  attached  to  circuits 
which  were  mounted  on  a  2  foot  x  16  foot,  1/8  inch  thick  aluminum  ground  plane 
to  provide  the  terminal  loads  shown  in  Fig.  3—1.  The  wires  were  held  together 
with  nylon  lacing  cord.  Wire  25  was  driven  at  the  left  end  by  a  sinusoidal 
oscillator  whose  frequency  was  varied  from  1  kHz  to  100  MHz  in  discrete  steps. 
The  other  end  of  the  wire  was  terminated  to  the  ground  plane  with  a  resistor  R. 
Wire  14  was  also  terminated  at  both  ends  with  a  resistor  R.  The  received 
voltage  (crosstalk)  at  the  left  end  of  this  wire  was  measured.  The  voltage 
transfer  ratio  is  then  the  ratio  of  this  received  voltage  to  the  voltage  ap¬ 
plied  to  wire  25.  Two  values  of  R  were  used:  R  =  50  ft  and  R  =  1  kft„  Those 
values  were  chosen  so  that  the  two  coupling  mechanisms  -  capacitive  and  induc¬ 
tive  coupling  -  could  be  accentuated.  The  characteristic  impedance,  Z  ,  of 

v 

one  #2 2  gauge  wire  2.4  cm  above  ground  is  on  the  order  of  300  0.  Thus  when 
R  =  50  ti,  R  <  Z„  and  capacitive  coupling  will  predominate.  When  R  =  1  kO, 

R  >  Z  and  inductive  coupling  will  predominate. 

In  order  to  investigate  the  sensitivity  of  the  crosstalk  to  variations  in 
relative  wire  position  as  well  as  load  impedance,  the  following  experiment  was 
conducted.  The  bundle  was  initially  wrapped  and  suspended  2.4  cm  above  the 
ground  plane  by  styrofoam  supports.  Next,  the  R  =  50  0  loads  were  attached 
and  the  crosstalk  measured  for  frequencies  from  1  kHz  to  100  MHz.  Then  the 
R  =  1  kfi  loads  were  attached  and  the  frequency  response  remeasured.  These  data 
constitute  the  results  for  cable  #1. 
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The  lacing  cord  was  then  removed,  the  cable  rewrapped  and  the  above  data 


were  taken  again.  This  constitutes  the  data  for  cable  if 2.  This  was  repeated 
8  times  to  give  data  for  a  total  of  10  cables.  Note  that  the  R  =  50  S2  and 
R  =  1  kfl  data  for  each  cable  were  for  the  exact  same  cable.  Thus  any  dif¬ 
ferences  in  sensitivity  which  are  found  for  these  two  different  loads  would 
clearly  indicate  the  dependence  on  load  impedance  of  that  sensitivity. 

The  results  are  plotted  in  Fig.  3-3  for  R  =  50  "1  and  in  Fig.  3-4  for 
R  =  1  kfi.  Note  in  Fig.  3-3  that  for  R  =  50  Q,  there  is  virtually  no  sensitiv¬ 
ity  to  variations  in  relative  wire  position  caused  by  rewrapping  the  cable 
bundle.  On  the  other  hand,  note  in  Fig.  3-4  for  R  =  1  kft  that  a  change  in 
crosstalk  of  as  much  as  20  dB  occurs  when  the  cable  bundle  is  rewrapped. 

Note  also  that  this  occurs  even  for  frequencies  where  the  cable  is  electrically 
very  short.  Thus  for  R  =  1  kfi,  these  10,  supposedly  identical  cable  bundles 
in  fact  exhibit  quite  different  crosstalk  characteristics. 

Also  plotted  on  these  data  are  the  predictions  of  a  typical  computer 
program,  the  Intrasystem  Electromagnetic  Compatibility  Analysis  Program, 

IEMCAP  [21].  This  is  representative  of  similar  "worst  case"  models  which 
attempt  to  bound  this  variability.  In  addition,  upper  and  lower  bounds  ob¬ 
tained  with  tolerance  intervals  based  on  the  statistics  of  the  data  are  shown 
[20]. 

The  crosstalk  data  may  represent  an  important  part  of  some  overall  trans¬ 
fer  function  for  interference.  Even  though  the  remainder  of  the  transfer 
function  is  characterized  by  little  or  no  variability,  variability  of  the 
crosstalk  may  render  any  deterministic  predictions  meaningless  unless  they  are 
upper  bounds.  Even  then,  these  upper  bounds  may  be  severely  pessimistic  re¬ 
sulting  in  considerable  overdesign.  A  more  appropriate  method  of  characteriz- 


Crosstalk  (dB)  vs.  Frequency 


Fig.  3-3.  Data  for  R  =  50Q  loads. 


ing  this  crosstalk  for  the  purposes  of  estimating  potential  problems  would 
seem  to  be  the  statistical  approach  alluded  to  in  the  previous  chapter. 


Other  portions  of  an  overall  transfer  function  may  exhibit  similar 
variabilities.  A  specific  example  of  this  is  discussed  in  the  next  chapter 
with  regard  to  variability  in  the  nonlinear  transfer  functions  of  Op  Amps. 

In  chapter  5,  an  experiment  is  detailed  which  contains  both  aspects  of  this 
variability.  The  purpose  will  be  to  investigate  the  feasibility  of  combining 
separate  sets  of  statistics  on  various  portions  of  a  transfer  function  to  de¬ 
termine  the  overall  statistics  of  the  transfer  function. 
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IV.  Statistics  of  Measured  Demodulation  RFI  responses 

in  Operational  Amplifier  Circuits. 

The  specific  RFI  effect  investigated  was  demodulation  of  amplitude 
modulated  (AM)  RF  signals  in  operational  amplifiers  (op  amps) 
to  produce  undesired  low  frequency  responses  at  the  AM-modulat ion 
frequency.  The  undesired  demodulated  response  may  then  be  processed 
as  a  desired  low  frequency  signal  by  the  low  frequency  compon¬ 
ents  that  follow  the  op  amp.  In  this  chapter,  the  investigation  of 
an  op  amp  unity  gain  buffer  amplifier  configuration  and  a  3-stage  [2] 
op  amp  configuration  will  be  discussed.  The  op  amp  types  chosen  are 
the  widely  used  741  bipolar  op  amps  which  have  conventional  npn  input 
transistors,  the  newer  LM10  bipolar  op  amps  which  have  less  conven¬ 
tional  pnp  input  transistors,  the  LF355  JFET-bipolar  op  amps  which 
have  junction  field-effect-transistor  (JFET)  input  transistors,  and 
the  CA081  MOS-bipolar  op  amps  which  have  metal-oxide-semiconductor 
field-effect  transistor  (MOSFET)  input  transistors.  Information  on 
the  units  tested  is  summarized  in  Table  4-1  and  Table  4-2.  The  second- 
order  transfer  function  l^C-f  ,f  ),  which  characterizes  the  RFI  demod¬ 
ulation  response,  was  determined  from  the  measured  data.  The  scatter 
plot  and  statistics  such  as  mean  and  standard  deviation  of  ^(-f  ,  f  ) 


for  25  to  30  units  of  each  op  amp  type  tested  will  be  presented. 


TABLE  4-1 


OPERATIONAL  AMPLIFIERS  TESTED 


Desig. 

741 

LM10 

CA081 

LF355 

Type 

Bipolar 

Bipolar 

Bi-MOS 

Bi-FET 

Input 

npn 

pnp 

MOSFET 

JFET 

Interior 

bipolar 

bipolar 

bipolar 

bipolar 

Units 

30 

25 

30 

30 

Codes 

5 

4 

2 

4 

Manu. 

2 

1 

1 

1 

TABLE  4-2 


MAKER  AND 

PRODUCTION 

DATA  INFORMATION 

OF  THE  30 

741  OP  AMPS 

Group 

Units 

Maker 

Year 

Week 

1 

5 

RCA 

1981 

27 

2 

5 

Fairchild 

1976 

36 

3 

10 

Fairchild3 

1981 

38 

4 

5 

Fairchild3 

1981 

35 

5 

5 

Fairchild'3 

1982 

49 

Manufactured  In  Hong  Kong. 
Manufactured  in  Korea. 


The  linear  response  of  an  amplifier  is  also  called  its  first- 


order  transfer  function.  Before  starting  measurements  of  the 
second-order  transfer  function  of  a  nonlinear  circuit,  measurement 
of  the  linear  response  is  usually  performed  to  check  the  normal 
operation  of  the  amplifier  circuit.  Figure  4-1  shows  the  experi¬ 
mental  set-up  for  measuring  the  linear  response  of  an  op  amp  unity 
gain  buffer.  In  this  measurement,  the  input  voltage  amplitude 
was  held  constant  at  50  mV  while  the  output  voltages  Vq^  and 
were  measured  corresponding  to  Rx  =  0  and  Rx  =  50  Q,  respectively. 
Then  the  linear  response  Vq/V^  corresponding  to  Rx  =  “  can  be 
derived  from  the  expression 


V 

O  ,100 

vi  = 


__50 _ ) 

V0l(mV)' 


R 

X 


co 


(4-1) 


The  reason  for  using  this  scheme  is  that  the  linear  response 
of  an  op  amp  unity  gain  buffer  is  a  function  of  load  resistance 
at  RF  frequencies  where  the  loop-gain  is  small.  The  coaxial  cable 
which  connects  the  output  node  Vq  to  the  input  of  RF  voltmeter  has 
to  be  terminated  by  a  50  Q  resistor,  which  loads  the  amplifier  out 
put.  Therefore  the  linear  response  corresponding  to  R^  =  ®  which 
corresponds  to  not  connecting  the  RF  voltmeter,  has  to  be  derived 
indirectly  from  Eq.  (4-1).  Figures  4-2  to  4-5  show  the  linear 
responses  of  the  unity  gain  configuration  of  the  741,  LM10,  CA081, 
and  LF355  op  amps,  respectively. 


Lineax  response  V  /v.  vs  frequency  for  U-A741  op  amp  #56  in  unity  gain  buffer  amplifier  circuit 


response  V  /V  vs  frequency  for  CA081  op  amp  #19  in  unity  gai 


Linear  response  V  /V,  vs  frequency  for  LF355  op  amp  #33  in  unity  gai 


Examining  the  data  shown  in  Figures  4-2  to  4-5,  we  observe 
that  the  unity  gain  buffer  stage  linear  frequency  response  is  down 
3  dB  at  0.03  MHz  for  the  LM10,  is  down  3  dB  at  0.7  MHz  for  the  741, 
is  down  3  dB  at  4  MHz  for  the  CA081,  and  is  down  3  dB  at  7  MHz  for 
the  LF355.  These  results  indicate  that  the  gain-bandwidth  product 
for  the  unity  gain  buffer  amplifier  configuration  is  lowest  for  the 
LM10  and  highest  for  the  LF355.  However,  we  should  note  that  the 
linear  responses  were  measured  for  only  one  unit  of  each  op  amp 
type.  Our  results  are  consistent  with  the  typical  data  for  these 
op  amps  provided  by  manufacturers  except  that  for  the  order  of  the 
CA081  and  LF355.  See  Ref.  [2]  for  a  plot  of  manufacturer's  typical 
data. 

4.2  Measurement  of  Demodulation  RFI  Response  of  Unity  Gain  Buffer 

From  the  previous  section,  we  have  just  shown  that  all  four 
op  amps  are  best  operated  at  frequencies  less  than  10  MHz  in  the 
unity  gain  buffer  configuration.  Now,  we  are  going  to  demonstrate 
how  amplitude-modulated  (AM)  RF  signals  are  demodulated  to  pro¬ 
duce  undesired  audio  frequency  signals  which  fall  into  the  normal 
frequency  range  of  the  op  amp  circuit  [22]. 

The  block  diagram  of  the  experimental  system  for  measuring 

the  demodulation  RFI  response  is  shown  in  Figure  4-6.  The  RF 

signal  generator  injects  into  the  non-inverting  input  of  the  op 

amp  circuit  an  amplitude-modulated  RF  signal,  V  ,  which  may  be 

8 

expressed  as 


of  the  RF  signal  at  the  receiving  end  of  the  cable. 


Experimental  set-up  for  demodulation  RFI  response  measurement  (after  Fang 
resistor  is  often  replaced  by  a  short-circuit  in  voltage  follower  circuits 


The  voltage  reading  at  the  tuned  AF  voltmeter  is  a  direct 
measure  of  the  demodulation  RFI  response  which  can  be  character¬ 
ized  by  the  second-order-transfer  function  H2(-fc>fc)  t*ie 
circuit  under  consideration.  An  expression  from  which  ^(-f^f^ 
values  can  be  determined  from  the  tuned  AF  voltmeter  readings 
has  been  derived  in  Chapter  2.  Using  Eq.  (2.3-36),  the  amplitude 
of  the  intermodulation  signal  at  frequency  f  is  given  by 

vm=mA2  |H2(-fc,fc)|  (4-3) 

The  AF  voltmeter  indicates  the  rms  voltage  of  the  AF  signal. 
If  we  denote  this  rms  voltage  by  VM,  then 

VM  =  0.707mA2  | H2(-fc,fc)j  (4-4) 

Equation  (4-4)  can  be  expressed  in  dB  with  respect  to  a  1  mV 
reference  level. 

20-log1()|VM/lmV|  -  20.1og10[0.707mA2|H2(-fc,fc)|]  -  20*log1()(10'3) 
=  57  +  20*log1Q(m)  +  40*logl0(A) 

+  20.1og10|H2(-£c,fc)l 


(4-5) 


The  amplitude  of  the  RF  carrier  voltage  A  can  be  given  in 


terms  of  the  generator  available  power  P  which  by  definition 


is  the  power  the  generator  would  deliver  to  a  load  =  Rg  where 

R  is  the  generator  internal  impedance.  Using  the  equivalent  cir- 
S 

cuit  of  Figure  4~8  with  R^  =  Rg,  we  obtain  the  relationship 


VT  2 

P  =  ~  =  — — 

gen  R±  8Rg 


(4-6) 


Note  that  is  an  rms  voltage,  but  that  A  is  an  amplitude.  If 


R^  =  Rg  =  50  ft,  Equation  (4-5)  can  be  expressed  in  terms  of  dBm 


as 


Pgen(dBm)  =  10*log10((A2/400)/10  3)  =  4  +  20*log1()(A)  (4-7) 

With  m  =  0.5,  Equation  (4-4)  can  be  written  as 

20.1og1()|VM/lmV|  =  2Pgen(dBm)  +  43  +  20-log^  |H2  (-fc,fc)  I  (4-8) 

Thus,  we  have  an  expression  relating  ^(-f^f  )  values  and 

the  tuned  AF  voltmeter  readings.  Besides,  Equation  (4-8)  also 

suggests  a  convenient  way  of  checking  whether  the  AF  voltmeter 

reading  are  caused  by  the  second-order  nonlinearities  of  the 

amplifier  circuit  under  test.  For  every  1  dBm  variation  in  Pgen> 

a  corresponding  2  dB  variation  should  be  observed  in  the  measured 

V„  value.  We  assume  this  condition  when  we  use  Equation  (4-8)  to 
M 

calculate  H„(-f  ,f  ).  Figure  4-9  shows  values  of  the  measured  AF 
Zee 

voltmeter  reading  in  dBV  versus  the  generator  available  power 
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Signal  Generator 


Amplifier  Under 
Test  with  Load 


Fig.  4-8.  Circuit  used  to  relate  the  available  power  P  that 

gen 

the  signal  generator  can  deliver  to  a  matched  load 
R.  =  R  to  the  signal  generator  voltage  amplitude  V  . 

O  © 


Pgen  for  four  types  of  op  amps  at  f£  =  10  MHz.  It  is  seen  that 
for  Pgen  <  0  dBm,  all  four  op  amps  have  data  plots  that  are 
straight  lines  with  slopes  equal  to  two  in  agreement  with  Equa¬ 
tion  (4-8).  The  region  in  which  this  behavior  is  observed  is 
often  called  the  square-law  response  region.  For  Pgen  >  0  dBm, 
the  slopes  of  the  data  curves  start  to  deviate  from  two  because 
of  nonlinear  terms  of  order  higher  than  two.  Whereas  we  cannot 
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Fig.  4-9.  Measured  values  of  AF  voltage  VM(rms)  at  the  tuned  AF 


voltmeter  vs  the  generator  available  power  P  for  the 


unity  gain  buffer  circuit  for  four  types  of  op  amps  at 
one  RF  frequency  : 


f  =  lkHZ ,  f  =  10MHz ,  m  =  0.5 
m  c 


set  Pgen  to  values  too  large,  it  is  also  not  appropriate  to  set 
it  to  values  too  low.  We  must  not  let  the  AF  signal  produced  by 


second-order  nonlinearities  of  the  amplifier  fall  below  the 

ambient  noise  level  and  become  difficult  to  measure.  Each  curve 

shown  in  Figure  4-9  is  used  to  obtain  a  H„(-f  ,f  )  value  at  one 

2  c  c 

RF  frequency.  Values  of  H2^“^c’^c^  at  ot^er  RF  frequencies  may 
be  obtained  from  plots  similar  to  that  shown  in  Figure  4-9.  It 
is  not  necessary  to  measure  as  many  data  points  as  shown  in 
Figure  4-9.  Often,  measurements  were  made  at  two  values  of  P 

gen 

such  as  -10  dBm  and  -20  dBm  to  verify  the  square-law  response. 


Then  the  value  of  VM(rms)  at  P  *  -10  dBm  was  used  to  determine 

M  gen 


H2('fc’fc 


). 


Statistical  Results  for  Unity  Gain  Buffer 


As  discussed  in  the  previous  section,  the  demodulation  RFI 
in  the  square  law  response  region  can  be  characterized  by  the  one 
parameter  HjC-f  ,f  ).  Shown  in  Figures  4-10  to  4-13  are  values  of 
^(-f  ,  fc)  for  RF  frequencies  ranging  from  0.1  to  400  MHz  for  the 
115  op  amps  tested.  In  Figure  4-10  data  for  30  741  op  amps  are 
presented.  Among  these  741  op  amps,  five  were  made  by  RCA  in  1981 
and  the  others  were  made  by  Fairchild  during  1976  to  1982.  The 
detailed  information  is  listed  in  Table  4-2.  The  H0(-f  ,f  )  values 
at  a  specific  RF  frequency  vary  +  1.5  to  +  11.5  dB.  In  Figure  4-11 


SECOND  ORDER  TRANSFER  FUNCTION  H2  VS  FREQUENCY  f 


is  of  the  second-order  transfe: 
frequency  for  thirty  CA081  (B: 


circuit  vs  RF  f 


data  for  25  LM10  op  amps  are  presented.  All  25  op  amps  were 
made  during  1981  and  1982  by  National  Semiconductor.  The  most 
interesting  aspect  of  the  LM10  data  is  that  two  LM10  units 
(called  mavericks)  have  ^(-f  , f£)  va^-ues  verY  different  from  the 
other  LM10  units.  It  should  be  noted  that  one  maverick  LM10 
op  amp  has  a  normal  linear  response  as  shown  in  Figure  4-18  of  Ref. 
[2].  Its  voltage  gain  vs  frequency  response  appeared  to  be 
similar  to  the  other  23  LM10  op  amps.  Since  the  maverick  LM10  might 
perform  as  expected  in  some  circuits,  its  data  can  not  really  be 
excluded  from  our  statistical  calculations.  The  23  similar  LM10  op 
amps  have  H„(-f  ,  f  )  values  varying  +  3  to  +  12  dB  at  any  one  RF  fre- 
quency.  In  Figure  4-12  data  for  30  CA081  op  amps  are  presented. 

All  30  units  were  made  by  RCA  during  1980.  The  ^(-f^f  )  values 
vary  +  1  to  +  5  dB  at  any  one  RF  frequency.  In  Figure  4-13  data 
for  30  LF355  op  amps  are  presented.  All  30  units  were  made  by 
National  Semiconductor  during  1980  to  1982.  The  most  interesting 
aspects  of  the  LF355  data  are  the  resonances  near  4.5  MHz  and  12 
MHz  where  the  H2(-fc>fc)  values  decrease  as  much  as  25  dB  and  15 
dB,  respectively.  The  cause  of  the  resonances  has  been  investi¬ 
gated  by  computer-aided  analysis  of  the  circuit  with  the  LF355 
op  amp  replaced  by  its  macromodel  as  described  in  Ref.  [2],  The 
H2(-fc,fc)  values  outside  the  resonance  regions  vary  from  +  3  dB 
to  +  8  dB  at  any  one  RF  frequency. 


Shown  in  Figures  4-14  and  4-15  are  the  mean  values  for 

which  are  denoted  by  H2  and  the  standard  deviation  a 
for  the  four  types  of  op  amps  tested.  Values  for  and  a  were 
calculated  using  the  data  for  all  the  op  amps  given  in  Figures 
4-10  to  4-13. 

The  mean  values  for  HgC-f^f  )  indicate  clearly  that  de¬ 
modulation  RFI  effects  are  greater  in  op  amps  with  bipolar  input 
transistors  (741  and  LM10)  than  they  are  in  op  amps  with  FET  in¬ 
put  transistors  (CA081  and  LF355).  This  is  a  most  important 
observation.  At  RF  frequencies  above  10  MHz,  demodulation  RFI  effects 
in  the  741  op  amp  are  significantly  greater  than  in  the  IH10  op  amp. 

It  is  believed  that  this  is  a  result  of  the  cutoff  frequency  of  the 
npn  bipolar  input  transistors  in  the  741  op  amp  being  higher  than 
the  cutoff  frequency  of  the  pnp  bipolar  input  transistors  in  the 
LM10  op  amp.  It  is  also  noted  that  the  mean  value  H2  for  the 
LF355  op  amp  shows  much  less  resonant  excursion  than  do  individual 
units.  This  is  a  result  of  the  spread  in  values  for  the  resonant 
frequencies  which  causes  a  smoothing  effect  for  For  RF  fre¬ 

quencies  less  than  10  MHz,  the  mean  value  is  lowest  for  the 
CA081  op  amps  (MOSFET  input  transistors).  For  RF  frequencies 
greater  than  10  MHz,  the  mean  value  ^  is  lowest  for  the  LF355  op 
amps  (JFET  input  transistors). 


The  values  for  the  standard  deviation  o  vs  RF  frequency  are 
shown  in  Figure  4-15.  It  is  noted  that  for  the  CA081  op  amps, 
the  o  values  are  less  than  2  dB  except  at  0.1  MHz  where  H2  is  very 
small.  For  the  LM10  op  amps,  the  a  values  were  calculated  both 
for  all  op  amps  and  also  for  the  normal  op  amps  only. 

2  *5 

Standard  Deviation  =  o  =  £ (E (x  -  y)  )  *  n] 

where  n  =  total  number  of  items, 
x  =  item  value, 
y  =  mean  value  of  all  items. 

Omitting  the  ^(-f  , fc)  values  for  the  two  maverick  LM10  op  amps 
decreased  a  by  approximately  2  to  4  dB  below  2  MHz  and  by  approxi¬ 
mately  1  dB  or  less  above  2  MHz.  The  0  values  shown  in  Figure 
4-15  are  the  values  calculated  for  all  25  op  amps.  For  the  LF355 
op  amps,  the  a  values  are  large  especially  in  resonant  regions  near 
4.5  MHz  and  12  MHz.  The  0  values  are  also  large  above  200  MHz  where 
the  H„(-f  ,f  )  values  are  small. 


Fig.  4-15.  Experimental  standard  deviation 
unity  gain  buffer  circuit  vs  RF 
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We  have  presented  results  for  the  unity  voltage  gain  buffer 
amplifier  circuit  which  is  often  called  a  voltage  follower.  That 
circuit  configuration  is  noninverting  because  the  intended  signal 
is  injected  into  the  noninverting  input  denoted  by  (+) .  The  RFI 
signal  was  also  injected  into  the  noninverting  input.  In  this 
section,  we  report  the  results  of  an  investigation  of  the  invert¬ 
ing  op  amp  configuration  shown  in  Figure  4-16.  The  intended 
signal  voltage  gain  is  A^  =  Vq^/V-j-^  =  “R2/Rl  and  the  intended 
signal  input  impedance  is  Rl.  The  capacitance  C4  in  the  feedback 
path,  when  included,  provides  RFI  suppression  [23],  An  experi¬ 
mental  configuration  similar  to  that  shown  in  Figure  4-6  was  used 
to  measure  demodulation  RFI  for  RF  frequqncies  in  the  range  100  kHz 
to  400  MHz.  Measurement  results  will  be  presented  for  Rl  =  10  kft 
and  R2  =  100  kfi  for  C4  =  0  pF  and  C4  =  27  pF  for  35  units  of  741 
op  amps.  Thirty  units  are  identical  to  those  listed  in  Table  4-2; 
the  five  additional  741  op  amps  were  made  by  National  Semicon¬ 
ductor  in  1982. 


The  actual  circuit  used  in  the  measurement  of  demodulation 
RFI  effects  in  inverting  op  amp  circuits  is  shown  in  Figure  4-17. 
This  circuit  will  be  referred  to  as  the  3-stage  op  amp  LED  circuit 
hereafter.  The  3-stage  op  amp  LED  circuit  was  designed  to  corre- 
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Fig.  4-16.  Inverting  op  amp  circuit  with  voltage  gain  *  -R2/R1. 
The  capacitor  C4,  when  included,  provides  RFI 
suppression. 

spond  to  a  circuit  in  a  warning  light  system.  In  a  warning  light 
system,  a  tranducer  generates  a  signal  which  must  be  amplified  suf¬ 
ficiently  in  order  to  turn  on  a  warning  light  such  as  an  LED  (Light 
Emitting  Diode).  The  RFI  signal  is  injected  into  the  inverting 
input  of  the  1st  stage.  The  RFI  signal  is  an  AM-modulated  Rf  signal 
Because  of  nonlinearities  which  exist  in  analog  circuits  called 
linear  amplifiers,  a  demodulated  RFI  response  at  the  AM  modulation 
frequency  is  generated  within  the  1st  op  amp  stage.  The  demodulated 


RFI  signal  produced  in  the  1st  stage  is  amplified  by  the  2nd  stage 
and  by  the  third  stage  which  drives  the  Light-Emitting  Diode  . 

The  amplified  demodulated  RFI  signal  can  light  the  LED  to  produce  a 
visual  display  of  the  existence  of  RFI.  The  EMI  experiment  was  de¬ 
signed  to  demonstrate  how  two  sets  of  statistics  from  two  indepen¬ 
dent  experiments  can  be  combined  to  predict  the  statistics  of  a  com¬ 
bined  experiment.  One  experiment  (See  Chapter  3)  provided  statis¬ 
tics  on  cable  coupling  for  wires  terminated  in  resistive  loads.  The 
experiment  described  in  this  section  and  in  the  following  section 
provided  second-order  transfer  function  statistics  for  the  3-stage 
op  amp  LED  circuit.  The  goal  was  to  combine  the  two  sets  of  inde¬ 
pendent  statistics  and  to  compare  to  the  statistics  of  the  combined 
experiment.  The  second-order  transfer  function  of  the  inverting  op 
amp  configuration  of  Figure  4-16  can  be  related  to  the  second-order 
transfer  function  of  the  3-stage  op  amp  LED  circuit  by  accounting 
for  the  linear  voltage  gains  of  the  second  and  third  stage  and  the 
attentuation  of  the  circuit  between  the  RF  generator  and  input  of 
the  1st  stage. 

This  paragraph  will  describe  the  function  of  each  element  in 
the  circuit  of  Figure  4-17.  The  RF  signal  generator  produces  an 
RF  voltage  with  an  RF  carrier  frequency  f  in  the  range  0.1  to 
400  MHz.  The  RF  voltage  is  AM-modulated  with  modulation  index 
m  =  0.5  and  AM-modulation  frequency  f^  ■  1  kHz.  The  capacitors 
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RF  Voltmeter  hp  302A  Wave  AnaU 


C6,  C7,  and  Cl  and  resistors  RIO,  Rll,  and  the  50  ft  input  resis¬ 
tor  of  the  RF  voltmeter  form  a  high-pass  RC  filter  which  blocks 
dc  ,and  any  spurious  1  kHz  components  from  the  RF  generator.  The 
first  stage,  is  the  op  amp  under  test;  it  has  an  inverting  config¬ 
uration.  Its  second-order  nonlinearities  cause  an  undesired  de¬ 
modulation  of  the  amplitude-modulated  RF  signal  to  produce  an 
undesired  audio  modulation  frequency  component  at  the  output 
of  the  first  stage.  The  ratio  of  the  feedback  resistor  R2  to  the 
input  resistor  R1  provides  the  intended  linear  voltage  gain  A  ^  of 
the  first  stage  at  low  frequencies.  The  capacitor  C4,  when  included, 

serves  as  an  RFI  suppression  capacitor  123],  The  RF  voltmeter 

RF 

reads  the  rms  RF  voltage  which  is  one-half  the  rms  RF  voltage 
at  the  input  of  the  first  stage.  Note  that  a  50  ft  resistor  is 
used  to  terminate  the  50  ft  coaxial  cable  connecting  the  circuit  and 
the  RF  voltmeter.  Furthermore,  the  resistor  combination  of  RIO  in 
parallel  with  the  series  combination  Rll  +  the  50  ft  termination 
resistor  at  the  RF  voltmeter  input  provides  at  high  frequencies  a 
50  ft  terminating  resistance  to  the  50  ft  coaxial  cable  connecting 
the  RF  generator  and  the  circuit.  Between  the  first  and  second 
stages  there  is  a  low-pass  filter  formed  by  R3,  R4,  and  C2.  It 
attenuates  any  RF  components  from  the  first  stage,  but  it  allows 
the  audio  frequency  component  to  pass  through.  The  second  and 
third  stages  provide  linear  voltage  amplification  of  the  audio-fre- 
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quency  component  V^.  The  second  stage  is  an  inverting  amplifier 
with  a  voltage  gain  of  100  determined  by  the  ratio  R5/ (R3  +  R4) . 

It  is  critical  to  include  a  resistor  equal  in  value  to  the  feedback 
resistor  R2  from  the  noninverting  input  of  the  first  stage  to  ground. 
That  resistor  reduces  the  dc  offset  voltage  at  the  output  of  the 
first  stage.  The  first  stage  dc  offset  voltage  will  be  amplified 
in  the  second  stage  and  can  cause  saturation  at  the  output  of  the 
second  stage.  The  resistor  connecting  the  noninverting  input  of  the 
second  stage  to  ground  serves  a  similar  purpose,  but  it  is  less 
critical.  The  capacitor  C3  and  resistor  R7  form  a  high-pass  filter 
with  a  160  Hz  cutoff  frequency  to  block  any  remaining  dc  offsets 
at  the  output  of  the  second  stage  from  entering  the  3rd  stage. 

The  third  stage  is  a  noninverting  amplifier  with  a  voltage  gain 
of  10  determined  by  the  ratio  (R8  +  R9)/R8.  The  LED  goes  on  when 
the  rms  voltage  of  the  amplifier  demodulated  audio  frequency  com¬ 
ponents  at  the  output  of  the  3rd  stage  exceeds  a  threshold 
value.  The  threshold  value  for  a  purely  1  kHz  signal  was  typically 
1.1  V  (rms)  which  corresponds  to  a  peak  voltage  of  1.5  V.  There 
were  occasions  when  the  LED  went  on  and  the  1  kHz  voltage  component 
was  much  less  than  1.1  V  (rms).  In  these  cases  a  strong  2  kHz  volt¬ 
age  component  was  observed.  A  2  kHz  voltage  component  can  also  be  a 
second-order  demodulation  RFI  component.  It  results  from  a  beating 
or  mixing  (intermodulation)  of  the  upper  and  lower  sidebands  of  the 


AM-modulated  RF  input  signal  in  the  1st  op  amp  stage.  Usually,  the 
2  kHz  RFI  component  was  in  the  range  0.15  +  0.005  V  (rms)  when  the 
1  kHz  RFI  component  was  1.0  V  (rms).  However,  when  a  resonant  effect 
was  observed  which  caused  the  1  kHz  voltage  component  to  be  suppressed, 
the  RF  voltage  has  to  be  increased  to  an  unusually  high  level 
in  order  to  turn  on  the  LED.  Then  a  strong  2  kHz  voltage  component 
was  observed.  The  2  kHz  voltage  component  might  be  as  large  as  the 
1  kHz  voltage  component.  Occasionally,  the  2  kHz  voltage  component 
exceeded  the  1  kHz  voltage  component.  We  suspect  that  these  circum¬ 
stances  indicated  that  a  4th-order  nonlinear  effect  may  have  mani¬ 
fested  itself. 

A  frequency  counter  was  used  to  measure  accurately  the  RF  fre¬ 
quency.  A  eathode  ray  oscilloscope  (CR0)  was  used  to  monitor  the  dc 
offset  voltage  and  the  AF  voltage  waveform  at  the  output  of  the  second  stage 

Prior  to  the  measurement  of  demodulation  RFI  responses,  the 
linear  voltage  gain  at  audio  frequency  1  kHz  of  the  entire  3-stage 
op  amp  LED  assembly  and  of  each  section  was  measured  to  verify  the 
circuit  connections.  An  Audio  Frequency  (AF)  oscillator  (hp651B) 
was  connected  in  place  of  the  RF  voltmeter  to  apply  an  AF  signal 
to  the  circuit  while  the  RF  generator  was  turned  off.  The  AF  volt¬ 
ages  V^,  V^2»  Vq^»  an<*  Vi^  w^ich  are  shown  in  Figure 

4-17  were  read  with  a  tuned  AF  wave  analyzer  (hp302A).  The  super- 


script  AF  denotes  Audio  Frequency  component.  The  AF  voltage 
at  the  output  of  the  3rd  stage  was  set  at  1.0  V  (rms)  while  readings 
of  other  AF  voltages  were  taken.  The  results  are  given  in  Table 
4-3.  The  reading  drew  our  attention.  When  the  op  amp  used  in 
the  2nd  stage  was  a  741  op  amp,  the  voltage  was  0.25  mV  instead 
of  the  expected  0.2  mV.  When  the  741  op  amp  was  replaced  with  a 
LF355  op  amp,  the  voltage  was  0.2  mV  in  agreement  with  what  we 
expected.  We  believe  that  the  discrepancy  in  the  V^  value  when  a 
741  op  amp  was  used  in  the  second  stage  is  related  to  the  low  input 
impedance  of  the  741  op  amp.  Other  AF  voltages  in  Table  4-3  have 
values  in  satisfactory  agreement  with  calculated  values. 


4.5  Measurement  of  Demodulation  RFI  Response  of  3-Stage  Op  Amp 

LED  Circuit 

In  all  measurements  of  demodulation  RFI  responses,  the  op  amps 

for  the  2nd  and  3rd  stage  of  the  3-stage  op  amp  LED  circuit  were 

the  same  (741-//20  and  741-//12  respectively).  The  1st  stage  op  amp 
was  changed.  Each  of  the  30  units  of  741  op  amps  listed  in  Table 
4-2  and  the  five  additional  741  op  amps  made  by  National  Semiconductor 

was  used  in  the  1st  stage.  The  signal  from  the  RF  generator  was  a 

50%  amplitude-modulated  (1  kHz)  RF  signal  as  before.  The  input  RF 
signal  at  a  specific  RF  frequency  was  adjusted  so  that  the  demodulated 
1  kHz  AF  component  at  output  V3  of  the  3rd  stage  reached  rms  values 
of  V^  =  1.0  V  and  0.4  V  sequentially.  At  the  same  time,  the  meter 


readings  of  the  RF  generator  and  RF  voltmeter  were  recorded.  This 


TABLE  4-3 


LINEAR  AF  NODE  VOLTAGES  AT  1  KHZ  OF  THE  3- STAGE  OP  AMP  LED  CIRCUIT 


hp  302A 

Expect) 

AF  Voltage 

Reading3 

Value 

^3  <V) 

1.0 

1.0 

^2  ("V) 

102 

100 

v£  <"v> 

0. 25b 

0.20 

C  <»s> 

0.98 

1.0 

vf  (yv) 

98 

100 

vf  (yv) 

115 

>  100 

All  readings  are  in  rms  values.  The  741  op  amps  were 
used  in  all  3  stages  (1st:  741-//1,  2nd:  741-//20, 

3rd:  741-//12). 

k  When  a  LF355  was  used  in  2nd  stage,  the  reading  was  0.2  mV. 

procedure  was  used  to  verify  the  square  law  response  region  described 
in  Section  4.2.  To  understand  this,  note  that  the  two  AF  voltages 
(1.0  V  and  0.4  V)  at  output  V3  corresponds  to  a  ratio  of  -8  dB. 
Therefore,  the  two  corresponding  meter  readings  of  the  RF  voltmeter 
or  RF  generator  should  give  a  ratio  of  approximately  -4  dB.  If 
otherwise,  the  data  taken  do  not  reflect  the  characteristics  of 
second-order  nonlinearities  and  should  be  rejected  for  use  in  deter¬ 
mining  H^.  If  the  data  correspond  to  the  square  law  response  region. 


the  second-order  transfer  function  can  be  determined  from  the  data. 


The  1  kHz  AF  rms  voltage  V^F  at  the  output  is  related  to  the 
RF 

RF  rms  voltage  V  from  the  RF  generator  by  the  second-order  transfer 
S 

function  H2(-fc,fc). 

/TV^  =  m(/Tv^)2  *lH2(-fc,fc)|  (4-9) 

where  m  is  the  modulation  index.  The  factor  2  arises  because 
amplitudes  are  used  in  nonlinear  transfer  function  expressions. 

From  Figure  4-17,  we  note  the  voltage  relationship  given  by 

VRF  =  4VRF  (4-10) 

g  m 

Substitute  Eq.  (4-9)  into  Eq.  (4-10)  and  solve  for  H  (-f  ,f  )  with 

M  C  C 

m  =  0.5.  The  result  is 

l“2(-fc’fc)l  -  ^3/(4Vf)2  W*U> 

AF  RF 

where  V__  and  V  are  in  rms  volts.  Equation  (4-11)  can  be 
01  tn 

expressed  in  dB  as 

Af  Rf 

201og10|H2(-fc,fc)|  -  201og10V03  -  401»g10Vn  -  21.0  (4-12) 

» 

The  second-order  transfer  function  ^  of  the  1st  stage  alone,  i.e., 
between  the  output  node  of  the  1st  stage  and  the  node  connecting  Cl 
and  Rl,  can  be  related  to  the  second-order  transfer  function  H2  of 
the  complete  circuit.  Accounting  for  the  linear  gains  of  the  2nd 
stage  (100)  and  the  3rd  stage  (10)  and  the  attenuation  factor  (0.5) 
of  the  input  attenuator,  we  obtain  the  results 
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lH2(“fc*fc) I  “  |H2(-fc,fc)|(2)Z/(100-10)  -  |H2<-fc,fc) |/250  (4-13) 

or  201og10|H2(-fc,fc)|  -  201og10|H2(-fc,fc)|  -  48.0  (4-14) 

4.6  Statistical  Results  for  3-Stage  Op  Amp  LED  Circuit 

As  discussed  previously,  the  demodulation  RFI  at  1  IcHz  in  the 
square  law  response  region  can  be  characterized  by  the  one  parameter 
H2(-fc,fc).  Shown  in  Figures  4-18  and  4-19  are  measured  values  of 
H2(-fc,fc)  of  the  3-stage  op  amp  LED  circuit  for  RF  frequencies  in 
the  range  of  0.1  to  150  MHz  for  35  units  of  741  op  amps  tested.  In 
Figure  4-18,  the  resistor  and  capacitor  combinations  are  R1  *  10  kf), 

R2  »  100  kft,  and  C4  =  C5  =  0.  The  H2(-fc,fc)  va^ues  at  a  specific 
RF  frequency  vary  +  3  to  +  11  dB.  In  Figure  4-19,  the  resistor  and 
capacitor  combinations  are  R1  =  10  kft,  R2  -  100  kfi,  and  C4  =  C5  =  27  pi 
The  H2^~^c’^c^  va*-ues  at  a  specific  RF  frequency  vary  +  4  to  +  15  dB. 
Whereas  H2  values  in  Figures  4-18  and  4-19  are  for  the  complete  3- 

f 

stage  op  amp  LED  circuit,  the  H2  values  for  the  1st  stage  alone  can 
be  obtained  easily  by  subtracting  48  dB  from  H2  values  for  the  com¬ 
plete  circuit.  See  Eq.  (4-14).  The  value  of  the  ordinate  in 

1 

Figures  4-18  and  4-19  is  reduced  by  48  dB.  The  variations  of  H2 
values  in  dB  at  a  specific  frequency  are  identical  to  those  of  H2 


values. 


741  Inverting  Amplifier 
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Shown  in  Figures  4-20  and  4-21  are  the  mean  values  for  H„(-f  ,f  ) 

C  C 

which  are  denoted  by  Hj  and  the  standard  deviation  a  for  35  units  of 
741  op  amps  tested  in  the  3-stage  op  amp  LED  circuit.  One  of  the 
two  sets  of  data  in  each  figure  corresponds  to  R1  ®  10  left,  R2  -  100 
kft,  and  C4  =  C5  =  0.  Another  set  corresponds  to  R1  =  10  kft,  R2  ■  100 

kft,  and  C4  =  C5  =  27  pF.  The  mean  values  for  ^(-f^f  )  indicate 

clearly  the  effect  of  the  RFI  suppression  capacitor  C4  when  it  is 
connected  in  the  feedback  path  of  the  1st  stage.  The  suppression 
of  demodulation  RFI  effects  caused  by  a  27  pF  capacitor  result  in  a 
reduction  in  from  3.5  dB  to  as  much  as  36.5  dB  at  a  specific  RF 
frequency.  Another  effect  of  the  RFI  suppression  capacitor  C4  is 
indicated  by  the  two  sets  of  standard  deviation  data  in  Figure  4-21. 
That  plot  shows  that  the  spreading  of  H2  values  is  wider  in  general 

by  including  C4  -  27  pF.  An  increase  as  much  as  5  dB  in  standard 

deviation  is  observed  at  certain  RF  frequencies. 

In  Figure  4-22,  RFI  demodulation  responses  in  two  741  op  amp 
circuits  are  shown.  One  set  of  data  are  mean  values  of  measured 
^  for  the  741  unity  gain  buffer.  Another  set  of  data  are  mean 
values  of  measured  ^  for  the  741  inverting  amplifier  with  a  voltage 
gain  of  10  and  no  RFI  suppression  capacitor  (C4  »  0).  For  compari¬ 
son  purposes,  the  values  for  the  741  unity  gain  buffer  from 

I 

Figure  4-14  have  to  be  adjusted.  We  want  an  for  the  unity  gain 

buffer  that  relates  the  audio-frequency  voltage  at  the  op  amp  out- 
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27  pF.  This  plot  shows  the  suppression  of  the  demodulation  RFI  effects 


i 


put  to  the  RF  voltage  at  the  noninverting  input  terminal.  We  must 

account  for  the  voltage  division  at  tne  output  which  contributes  a 

factor  of  (2)  and  the  voltage  division  at  the  input  which  contri- 

2 

butes  a  factor  of  (2)  .  See  Figure  4-7.  Therefore,  the  relation- 

I 

ships  between  H2  and  H2  are  given  by 

lH2(-fc‘fc)i  =  SlV-Wl  (4"15) 

201og10iH2(-fc’fc)!  =  201og10|H2(-fc,fc)|  +  18.1  (4-16) 

The  result  is  an  increase  of  18  dB  in  the  H2  values  plotted  in 
Figure  4-14.  Also  H2  values  for  the  741  inverting  amplifier  from 
Figure  4-20  have  been  adjusted  using  Eq.  (4-14).  The  result  is  a 
decrease  of  48  dB  in  the  values  plotted  in  Figure  4-20.  The  adjusted 

I 

H2  values  in  Figure  4-22,  therefore,  exclude  the  linear  amplification 
and  attenuation  stages  in  the  experimental  setups.  The  adjusted 

I 

mean  values  for  H2  are  surprisingly  similar  in  that  both  sets  of 
data  points  have  a  peak  value  near  10  dB.  The  comparison  shows  that 
the  741  unity  gain  buffer  has  a  lower  average  RFI  response  for  RF 
frequencies  below  8  MHz.  Above  8  MHz,  the  average  RFI  demodulation 
response  in  the  741  inverting  amplifier  is  lower. 
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In  the  previous  two  chapters,  examples  of  significant  variations  in  coup¬ 


ling  characteristics  of  a  transfer  function  were  shown.  The  purpose  of  this 
chapter  is  to  review  an  experiment  which  combined  these  two  transfer  functions, 
The  object  of  the  experiment  is  to  investigate  the  feasibility  and  adequacy  of 
combining  individual  statistics  on  separate  parts  of  the  overall  transfer  func¬ 
tion. 

5.1  Cable  Coupling  Variability 

In  Chapter  3,  the  sensitivity  of  cable  coupling  (crosstalk)  to  variations 
in  relative  wire  position  was  examined.  It  was  found  that  the  sensitivity  of 
crosstalk  to  variations  in  relative  wire  position  was  a  strong  function  of  the 
values  of  the  load  impedances.  It  was  found  that  for  "high  impedance"  loads, 

R  =  1  kft,  that  the  sensitivity  to  changes  in  wire  position  could  be  extraor¬ 
dinarily  large.  Changes  in  crosstalk  caused  by  rewrapping  the  25-wire  cable 
bundle  caused  variations  of  as  much  as  20  dB  even  for  frequencies  where  the 
line  was  electrically  very  short. 

The  25  wire  cable  bundle  described  in  Chapter  3  was  reconstructed.  Ad¬ 
ditional  data  were  taken  in  the  manner  described  in  chapter  3  for  R  =  50  ft  and 
R  =  1  kft  loads.  The  bundle  was  rewrapped,  loads  of  R  =  50  ft  were  attached  and 
the  crosstalk  was  measured  for  frequencies  of  4  MHz,  5  MHz,  10  MHz  and  15  MHz. 
Then  loads  of  R  =  1  kft  were  attached  and  the  crosstalk  data  retaken.  This 
constituted  the  data  for  bundle  11.  (The  10  data  points  obtained  previously 
constitute  cables  1  through  10.)  The  bundle  was  rewrapped  and  the  above  data 
retaken.  This  constitutes  cable  12. 

This  process  was  repeated  to  yield  data  for  a  total  of  60  cables.  The 
results  at  frequencies  of  4  MHz,  5  MHz,  10  MHz  and  15  MHz  are  shown  in  Fig. 


5.1-1  for  R  =  50  ft  and  for  R  =  1  kft.  Note  that  a  similar  sen¬ 
sitivity  uncovered  previously  is  found  for  these  additional  data.  For  R  =  50 
ft,  we  observe  virtually  no  sensitivity  to  changes  in  relative  wire  position 
caused  by  rewrapping  the  bundle.  There  is  more  sensitivity  at  15  MHz  than  at 
4  MHz,  5  MHz  and  10  MHz.  This  is  because  the  cable  is  becoming  electrically 
long.  At  15  MHz,  the  cable  is  approaching  1/4  X  in  electrical  length.  Thus 
sensitivity  to  variations  in  wire  position  can  be  expected  to  become  more 
severe  at  this  frequency  than  at  4  MHz,  5  MHz  and  10  MHz.  However,  note  that 
the  variability  of  the  crosstalk  caused  by  rewrapping  the  cable  bundle  with 
R  =  1  kft  is  significant  at  all  frequencies.  The  difference  in  the  maximum 
and  minimum  crosstalk  for  the  60  bundles  is  on  the  order  of  15  dB  at  4  MHz. 

These  data  will  be  used  in  a  combined  experiment  which  is  discussed  at 
the  end  of  this  chapter.  In  the  next  section  we  will  investigate  the  varia¬ 
bility  in  the  nonlinear  transfer  functions  of  Op  Amps. 
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Fig.  5.1-1.  Random  cable  sensitivity  data  for  R  =  50 ft  and  R  =  lkfi  loads. 
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5.3  The  Combined  Experiment 


The  data  of  the  previous  two  sections  have  shown  the  rather  extreme  vari¬ 


ability  which  can  occur  in  certain  transfer  functions  which  are  important  con¬ 


tributors  to  electromagnetic  interference.  In  this  section  we  will  combine 


those  variabilities  in  a  single  experiment.  Our  interest  will  be  to  obtain 


statistics  on  some  overall  transfer  function  which  has  these  two  individual 


parts  as  contributors. 


A  block  diagram  of  the  experiment  which  was  performed  was  shown  in  Chap¬ 


ter  Hand  is  repeated  in  Fig.  5.3-1.  The  25-wire  cable  bundle  described  in 


Chapter  3  is  driven  by  an  HP606-C  oscillator.  This  carrier  frequency  was  set 


at  4  MHz  and  15  MHz.  A  1  kHz  modulation  was  applied  to  this  carrier,  and  the 


level  of  modulation  was  set  at  50%.  This  composite  signal  is  applied  to  the 


right  end  of  wire  25.  The  signal  is  coupled  via  crosstalk  to  the  right  end  of 


wire  14,  A  printed  circuit  board  is  attached  to  the  end  of  this  wire.  The 


detailed  schematic  of  the  board  layout  is  shown  in  Fig.  5.3-2.  Photographs  of 


the  board  attached  to  the  end  of  the  cable  bundle  are  shown  in  Fig.  5.3-3. 


The  first  stage  is  an  inverting  stage  with  a  DC  gain  of 


—  =  10 
R1 


This  is  followed  by  a  low  pass  filter  consisting  of  and  C ^  with  a  3dB  cut¬ 


off  frequency  of  1764  Hz.  This  is  followed  by  another  inverting  amplifier 


having  a  DC  gain  of 


R3  +  R4 


This  second  stage  is  followed  by  a  highpass  filter  having  a  3dB  cutoff  frequen¬ 


cy  of  159  Hz.  A  noninverting  buffer  stage  with  a  DC  gain  of 
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Printed  Circuit  Board 


provides  the  final  level  to  the  LED.  Thus  the  overall  DC  gain  is 

V  V  V  V  V  V 
3  3  4  V5  6  7 

v~=\rx\rx\rxv“x\r 

2  v4  5  6  V7  z 

=  9898 

The  functioning  of  the  circuit  is  as  follows.  The  AM  modulated  signal  at 

is  coupled  via  crosstalk  to  the  input  of  the  board  appearing  at  V2.  After 

passage  through  the  first  stage,  DC,  1  kHz  and  2  kHz  components  appear  at  the 

output  of  that  stage  at  V^.  In  addition  any  AM  modulation  which  may  exist  is 

removed  by  the  low  pass  filter  resulting  in  a  baseband  signal  at  V  .  The 

6 

resulting  baseband  signal  is  passed  to  the  second  stage  which  provides  ampli¬ 
fication.  A  highpass  filter  eliminates  any  DC  which  may  be  due  to  nonlinear 
rectification  or  DC  offset  resulting  from  the  first  two  stages.  During  the 
course  of  the  experiment  the  DC  level  at  V,.  was  monitored  to  ensure  that  the 
second  stage  was  not  driven  into  saturation.  Finally,  the  resulting  baseband 
signal  at  is  passed  through  a  buffer  amplifier  (the  third  stage)  to  the  LED. 

A  total  of  52  741  Op  Amps  were  purchased.  Of  that  total,  24  were  National 
Semiconductor  LM741  and  28  were  Fairchild  pA  741.  Two  of  the  Fairchild  pA741 
units  were  used  in  the  second  and  third  stages  throughout  the  experiment  while 
the  remaining  50  units  were  used  in  the  first  stage. 

A  Hewlett-Packard  8405A  Vector  Voltmeter  was  used  to  measure  the  carrier 
level  at  V^  (channel  A).  (A  10:1  divider  was  used  on  the  A  channel  probe  to 
insure  that  it  would  not  be  damaged  by  voltages  greater  than  2V.)  A  voltage 
divider  was  constructed  on  the  board  at  V2  as  shown  in  Fig.  5.3-2.  This  con¬ 
sisted  of  a  series  combination  of  560  +  390  =  950ft  connecting  and  the  center 


pin  of  a  BNC  connector  mounted  on  the  board.  Channel  B  of  the  Vector  Volt¬ 
meter  was  connected  to  this  point  (V^)  through  a  50ft  feedthrough  resistor. 
Therefore  V^’  =  V2/2O.  (The  actual  ratio  is  19.8675.)  Thus  the  impedance 
from  V2  to  ground  was  approximately  1  kft.  (The  input  impedance  of  the  Vector 
Voltmeter  provided  negligible  loading  across  the  50ft  feedthrough.)  The  load 
impedance  seen  by  wire  14  at  the  input  to  the  board  was  then  1  kft  in  parallel 
with  the  impedance  seen  looking  into  the  first  stage.  The  total  input  impe¬ 
dance  looking  into  the  board  was  measured  with  an  HP  4815A  Vector  Impedance 
Meter.  At  4  MHz,  the  input  impedance  to  the  board  was  found  to  be  895/-12°  ft, 
whereas  at  15  MHz  this  impedance  was  found  to  be  720/-38°  ft.  The  output  of 
the  second  stage  at  V,.  was  monitored  with  a  Tektronix  2213  oscilloscope  to 
ensure  that  the  DC  levels  were  small  enough  so  that  A2  was  not  driven  into 
saturation.  In  all  cases  this  level  was  less  than  one  volt.  Voltage  V^  was 
monitored  with  an  HP  3400A  RMS  Voltmeter.  The  output  level,  V^,  was  monitored 
with  an  HP  3400A  RMS  Voltmeter  and  an  HP  3580A  Low-Frequency  Spectrum  Analyzer. 
The  reason  for  using  both  these  instruments  to  monitor  V^  was  that  V^  con¬ 
tained,  in  addition  to  1  kHz,  a  2  kHz  component.  The  RMS  Voltmeter  shows  the 
combination  of  these  two  components,  whereas  the  spectrum  analyzer  clearly 
separates  them.  It  was  found  that  in  all  cases  the  2  kHz  component  at  V^  was 
approximately  16  dB  down  from  the  1  kHz  component.  It  was  the  1  kHz  component 
at  V^  which  was  recorded. 

In  each  case  the  carrier  levels  at  V^  required  to  produce  1  kHz  levels  at 
Vj  of  .2V,  ,4V,  IV  and  sufficient  to  barely  light  the  LED  (normally  around  1.25 
V)  were  recorded.  The  procedure  is  as  follows: 

1)  Wrap  cable  bundle. 

2)  Insert  A1  Op  Amp  after  trimming  and  straightening 


3)  Turn  on  power  supply  for  Op  Amps. 

4)  Turn  up  modulation  to  50%,  turn  up  carrier  to  get 


V3  =  . 2V. 

HP3590A  set  on  linear  and  ,2V/div 

HP3400A  set  on  .3V/div 

HP3400A  V.  set  on  .03V/div 
4 

change  HP3580A  to  .lV/div  to  read  2  kHz. 

5)  Turn  down  modulation. 

6)  Read  V^/10  on  channel  A  of  vector  voltmeter. 

7)  Read  on  channel  B  of  vector  voltmeter. 

8)  Turn  up  modulation  to  50%  and  set  =  ,4V 
HP3580A  set  on  linear  and  lV/div 
HP3400A  V3  set  on  lV/div 

HP3400A  V4  set  on  .lV/div 

change  HP3580A  to  ,2V/div  to  read  2  kHz. 

9)  Turn  down  modulation. 

10)  Read  V^/10  on  channel  A  of  vector  voltmeter. 

11)  Read  on  channel  B  of  vector  voltmeter. 

12)  Turn  up  modulation  to  50%  and  set  V3  =  IV. 
HP3580A  set  on  linear  and  2V/div 
HP3400A  V3  set  on  lV/div 

HP3400A  set  on  .lV/div 

change  HP3580A  to  lV/div  to  read  2  kHz. 

13)  Turn  down  modulation. 


15)  Read  on  channel  B  of  vector  voltmeter. 

16)  Turn  up  modulation  to  50%  and  set  to  light 
(barely)  LED 

HP3580A  set  on  linear  and  2V/div 

HP3400A  V3  set  on  3V/div 

HP3400A  V4  set  on  .3V/div 

change  HP3580A  to  lV/div  to  read  2  kHz 

17)  Turn  down  modulation. 

18)  Read  V^/10  on  channel  A  of  vector  voltmeter. 

19)  Read  on  channel  B  of  vector  voltmeter. 

20)  Change  carrier  frequency  to  15  MHz  and  repeat 
steps  (1)  -  (19) 

21)  Rewrap  bundle  and  repeat  steps  (1)  -  (20). 

Thus  there  are  obtained  50  data  points  for  each  carrier  frequency  and 

level.  Each  data  point  gives  the  ratio  of  the  1  kHz  component  at  to 

the  (unmodulated)  carrier  level  at  for  one  bundle  wrap  and  one  of  the  50 

741  Op  Amps.  It  should  be  observed  that,  from  the  data  which  was  recorded, 

are  able  to  obtain  (1)  the  cable  coupling  transfer  function: 

(carrier) 

(carrier) 

as  well  as  (2)  the  nonlinear  transfer  function  of  the  board: 

V3  (1  kHz) 

V 2  (carrier) 

The  latter  transfer  function  is  related  to  ^  of  the  board. 

In  the  next  chapter,  these  data  will  be  analyzed. 


The  feasibility  experiment  consisted  of  three  independent  investigations. 
One  experiment  studied  the  crosstalk  variations  in  a  25-wire  cable  bundle 
caused  by  unwrapping  and  rewrapping  the  bundle.  A  second  experiment  examined 
fluctations  in  the  second-order  nonlinear  transfer  function  of  operational 
amplifiers.  The  third  experiment  investigated  variability  in  the  amplitude 
of  an  interfering  signal  needed  to  switch  on  an  LED  when  a  cable  bundle 
was  rewrapped  and  an  operational  amplifier  changed  each  time  the  experiment 
was  performed.  These  experiments  are  discussed  in  detail  in  Chapter  V.  In 
this  chapter  some  of  the  experimental  data  is  analyzed. 

6. 1  Analysis  of  the  Crosstalk  Data 

The  crosstalk  within  the  cable  bundle,  from  to  V2,  was  characterized 
by  the  linear  voltage  transfer  function,  G(f).  The  probability  density 
function  associated  with  | G (f ) [  was  denoted  by  f  (g)  .  The  data  obtained  for 
1 G(f ) |  is  listed  in  Fig.  5.1-1  as  a  function  of  load  resistance  and  frequency. 
Since  the  cable  input  voltage  was  1  volt,  the  cable  output  voltage  equals 
|G(f)|. 

For  analytical  convenience,  it  is  desirable  to  approximate  f  (g)  by  a 

u 

well  known  distribution.  In  general,  several  different  distributions  may 
provide  acceptable  approximations  to  the  experimental  data.  The  procedure  for 
finding  suitable  distributions  consisted  of  two  steps.  First,  the  coe ff icients 
of  skewness  and  kurtosis,  two  measures  which  contain  information  relative  to 
the  shape  of  a  distribution,  were  used  to  make  preliminary  selections  of  can¬ 
didate  distributions.  The  candidates  were  then  tested  to  see  whether  any 
should  be  rejected  on  the  basis  of  a  significant  statistical  deviation  between 


95 


the  experimental  data  and  the  distribution  being  tested.  Approximations  to 
fg(g)  were  then  chosen  from  those  distributions  which  passed  the  test.  How¬ 
ever,  it  was  not  possible  to  state  whether  one  distribution  which  passed  the 
test  was  a  better  approximation  than  another  distribution  which  also  passed 
the  test. 

til 

The  coefficients  of  skewness  and  kurtosis  are  related  to  the  kr-— 
central  moment  of  a  random  variable  where  k  is  a  positive  integer.  Let  the 
mean  of  the  random  variable  G  be  denoted  by  m  .  The  k— central  moment  of 
G  is  defined  to  be 


^Gk  = 


where  E[  ]  denotes  the  statistical  operation  of  expected  value. 


(6.1-1) 


definition. 


ou,  =  - TT'i  =  coefficient  of  skewness 

G3 


01  ,  *  - j  =  coefficient  of  kurtosis. 

[yG2J2 


(6.1-2) 


(6.1-3) 


The  coefficient  of  skewness  is  a  measure  of  the  asymmetry  of  the  probability 
density  function.  The  coefficient  of  kurtosis  is  a  measure  of  the  "peaked¬ 
ness"  or  "flatness"  of  the  probability  density  function  in  the  central  part  of 
the  distribution. 

Because  the  underlying  distribution  of  G  is  unknown,  it  is  necessary  to 
use  the  experimental  data  to  estimate  m^,  ^g3*  an<^  ^G4*  ^et  ^G  c*enote 

the  number  of  rewrappings  of  the  cable  bundle.  Also,  let  g^  denote  the 
measured  crosstalk  linear  transfer  function  for  the  k*"*1  rewrapping  of  the 


v- 


V-V-  .VI*'  . 

aV-V- 
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bundle.  Unbiased  estimates  for  mG,  and  are  then  given  by 


a°'  Ui  v 


(6.1-4) 


^G2  Ng-1  (gk  “  mG)  ’ 


(6.1-5) 


^G3  (Ng-1)(Ng-2)  (gk"mG)  ’ 

1  NG2  "g  ,  xs  .4  6NG-9  .2 

JG4  2  gk  mG  2  ^G2  ' 

(NG-1)(NG  -3NG+3)k=l  *  G  ng  -3Ng+3 


(6.1-6) 


(6.1-7) 


Using  Eqs.  (6.1-4)  through  (6.1-7)  in  Eqs.  (6.1-2)  and  (6.1-3),  estimates 
of  the  coefficients  of  skewness  and  kurtosis  become 


wG2r 


^G2-* 


(6.1-8) 


(6.1-9) 


Preliminary  selections  of  some  well  known  distributions  for  approximat¬ 
ing  f  (g)  are  made  by  comparing  the  estimates,  &  and  ,  with  the  known 
G  G3  G4 

values,  aG2  and  a^,  of  the  various  distributions  under  consideration. 

Those  distributions  for  which  a  reasonable  match  exists  become  candidates  for 
further  consideration.  In  general,  the  probability  density  functions  of  the 
well  known  distributions  contain  two  or  more  numerical  parameters  for  which 
numerical  values  must  be  determined  in  order  for  the  distribution  to  be  com¬ 
pletely  specified.  These  parameters  were  chosen  such  that  the  mean  and 

variance  of  the  distribution  equaled,  respectively,  the  sample  mean,  m  ,  and 

2 

the  sample  variance,  =  (1^.  can  shQWn»  however,  that  the 
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coefficients  of  skewness  and  kurtosis  are  independent  of  the  mean  and  vari¬ 
ance  of  a  distribution. 

Having  selected  a  group  of  possible  distributions  for  approximating 
fG(g)  on  the  basis  of  comparing  and  <5^  with  ando^,  respectively, 
it  is  then  desirable  to  perform  a  "goodness  of  fit"  test  on  each  candidate 
distribution.  This  enables  a  distribution  to  be  rejected  when  a  significant 
statistical  deviation  exists  between  the  experimental  data  and  the  distri¬ 
bution  being  tested.  Even  though  a  distribution  is  not  rejected  by  the  sta¬ 
tistical  test,  it  cannot  be  accepted  with  perfect  confidence.  In  addition, 
if  several  distributions  are  not  rejected,  the  goodness  of  fit  test  cannot 
be  used  to  accept  one  distribution  over  another.  The  best  that  can  be  said 
regarding  distributions  which  pass  the  test  is  that  they  were  not  rejected. 
The  goodness  of  fit  test  prevents  acceptance  of  a  distribution  which  is 

likely  to  be  a  poor  approximation  for  fG(g). 

2 

The  x  test  [16]  was  used  as  the  goodness  of  fit  test.  In  applying  this 
test,  the  range  of  G  was  divided  into  5  cells  such  that  it  was  equally  pro- 


The  significance  level,  3, 

o 


is  defined  to  be  the  conditional  probability  of 


rejecting  the  distribution  given  that  the  distribution  under  test  is  a 


suitable  approximation.  In  applying  the  x  test,  the  significance  level  was 


chosen  to  be 


In  addition, 


3  =  0.01 

O 


ng  -  60 


(6.1-12) 


(6.1-13) 


for  the  crosstalk  data  of  Fig.  5.1-1.  The  numerical  value  obtained  for  the 
2 

X  statistic,  using  Eq. (6.1-11),  was  then  compared  to  a  threshold  value  tabu¬ 
lated  for  Bo=  0.01  and  the  distribution  under  test  [16].  The  distribution 

under  test  was  accepted  as  a  suitable  approximation  to  f  (g)  provided  the 

G 

2 

computed  value  of  the  x  statistic  was  less  than  the  tabulated  threshold. 

The  data  of  Fig.  5.1-1  reveal  that  the  lkft  resistive  load  resulted  in 
significant  crosstalk  variability  at  all  frequencies  tested.  Consequently, 
in  the  feasibility  experiment,  the  25-wire  cable  bundle  was  terminated  in 
approximately  a  lkft  resistive  load.  The  carrier  frequency  of  the  AM  inter¬ 
ference  was  selected  to  be  4  MHz  and  15  MHz.  Therefore,  approximate  probabil¬ 
ity  density  functions  for  f „ (g)  were  found  only  for  the  data  corresponding  to 

G 

1)  R  =  lkfl,  fc  =  4  MHz 

2)  R  =  lkfl,  f  =  15  MHz. 

c 

This  data  is  found  in  the  first  and  fourth  columns  of  Fig.  5.1-1.  The  two 
cases  are  discussed  separately. 

Case  1:  R=  1  kT<i ,  f  =  4  MHz 

c 

Using  the  data  from  the  first  column  of  Fig.  5.1-1,  the  estimated  mean, 
variance,  coefficient  of  skewness,  and  coefficient  of  kurtosis  were  calculated1 
to  be 
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7.002  x  10  volts 


0 eQ )2  =  aG2  =  5.48  x  10-4  (volts)2 

aG3  «  0.859 

a  =  3.285 
G4 


(6.1-14) 


Since  the  coefficients  of  skewness  and  kurtosis  are  independent  of  mean  and 
variance,  normalized  distributions  with  zero  mean  and  unit  variance  were  con¬ 
sidered  as  candidates.  These  distributions  are  convenient  to  work  with  be¬ 
cause  only  one  parameter  is  needed  to  characterize  them  [16].  The  selected 
candidates  are  listed  below  along  with  their  coefficients  of  skewness  and 
kurtosis: 

a)  Log  Normal  (C  =  1/4),  =  0.778,  =  4.096 

b)  Weibull  (6  =  2),  =  0.63,  =  3.24 

c)  Gamma  (a  =  7),  =  0.71,  =  3.75. 

2 

Application  of  the  X  test  resulted  in  acceptance  of  the  Log  Normal  and 
Gamma  distributions  but  in  rejection  of  the  Weibull.  Inserting  the  sample 
mean  and  variance  from  Eqs . (6 . 1-14) ,  the  analytical  expressions  for  the  den¬ 
sity  functions  of  the  Log  Normal  and  Gamma  distributions  are  given  by: 


Log  Normal  (C  =  1/4) 


1 . 596 (g  -I-  .022)  1  exp{-8[£n((11.193)(g+.022))]2},  g_>  -.022 


fG(g) 


(6.1-15) 

g<  -.022 


Gamma  (a  =  7) 

f„(g) 


7, 


016  x  10)  12(g-.004)?  exp {  -120.831 (g-. 004) }  ,  g^  .004 


(6.1-16) 


The  two  density  functions  are  plotted  in  Figs.  6. 1-1 (a)  and  (b) ,  respectively 
Observe  that  multiplicative  factors  are  noted  along  the  axes. 


Case  2:  R  =  lkfl.  f  =  15  MHz 

- e - 

Using  the  data  fro,m  the  fourth  column  of  Fig.  5.1-1,  the  estimated  mean, 
variance,  coefficient  of  skewness,  and  coefficient  of  kurtosis  were  calculated 
to  be 


m  =  1.03  x  10  1  Volts 

(j 


aG2  =  4.30  X  10  3  (Volts)2 


SG3  "  K224 


aG4  =  4.279. 


(6.1-17) 


The  selected  candidates  are  listed  below  along  with  their  coefficients  of 
skewness  and  kurtosis: 

a)  Log  Normal  (C  =  1/4),  =  0.778,  =  4.096 

b)  Gamma  (a  =  3) ,  =  1,  =  4.5 

c)  Gamma  (a  =  4) ,  =  0.89,  =  4.2. 

2 

Application  of  the  X  test  resulted  in  acceptance  of  both  Gamma-distributions 
but  in  rejection  of  the  Log  Normal.  Inserting  the  sample  mean  and  variance 
from  Eqs.  (6.1-17),  the  analytical  expressions  for  the  density  functions  of 
the  Gamma  distributions  are  given  by: 


Gamma  (a  =  3) 

r 


fG^  ' 


(1.444  x  105) (g+.028)3  exp{-30.51  (g+.0281)}  ,  g  > 


.028 

(6.1-18) 

.028 


101 


/v;>7 

li-RiVy  A  .  -V*  .  -  , 


Gamma  (a  =  4) 


fG(8) 


4.924  x  106)(g  +  0.44)4  exp  {-34.111 (g+.044)  }  ,  g>-.044 


(6.1-19) 


»  g<~ .044 


The  two  density  functions  are  plotted  in  Figs.  6. 1-2 (a)  and  (b) , 


respectively . 


6.2.  Analysis  of  the  Second-Order  Nonlinear  Transfer  Function  Data. 


Demodulation  of  the  AM  interference  was  characterized  by  H2(-fc»fc), 


the  second-order  nonlinear  transfer  function  relating  the  intefering  signal 


at  to  the  demodulated  1  kHz  tone  at  .  In  determining  j (— f ^ , f j  > 

was  measured  directly.  However,  was  measured  using  a  2:1  voltage  divi¬ 
der.  If  the  measurements  for  and  are  denoted  by  and  V~m,  respec¬ 


tively,  then 


V2  "  2V2m 


(6.2-1) 


V3  "  V3m* 


From  Eq.  (2.4-2) 


A3  =  mA2  |H2(-fc,fc) [. 


(6.2-2) 


Therefore,  the  magnitude  of  the  second-order  nonlinear  transfer  function 


is  given  by 


(6.2-3) 


Recall  that  A^  and  A^  refer  to  peak  values.  Since  the  measurements  were 


recorded  as  rms  voltages, 


A2‘2^"  V2m 


A3  '  ^  V3» 


(6.2-4) 


*  v  v  v  v  v/  •  -w.v.v  v.  .  \  -\  «\  V*  .*  /•  /•  .\.s  S.-  Wv‘v  v 


(X10») 


(a) 


o 


(X 1 0  -») 


(b) 


Fig.  6.1-2.  Approximations  to  f  (g)  at  15  MHz.  (a)  Gamma  (a  *  3) , 
(b)  Gamma  (a  =  4) . 
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In  the  experiment  the  modulation  index  was  set  at 
m  =  0.5. 


(6.2-5) 


Consequently,  Eq. (6.2-3)  can  be  rewritten  as 


■Jl  V 


3m 


I H_(-f  » f  ) I  =  - 

2  (0.5)(2v/2  V2ffl) 


3m 


2  =  0.354  -^-2. 

V  2m' 


(6.2-6) 


Eq.  (6.2-6)  was  used  to  evaluate  [ H_ (— f  ,f  )|. 

L  C  C 


As  explained  in  section  5.2,  measurements  were  taken  on  both  pA  741 
and  LM  741  operational  amplifiers.  | H2 (— I  was  f°un<*  to  be  significantly 
smaller  at  15  MHz  for  the  LM  741  devices.  As  a  result,  it  was  not  possible 
to  light  the  LED  at  15  mHZ  when  LM  741’ s  were  inserted  into  the  first  stage 
of  the  PCB  during  the  feasibility  experiment.  Because  of  this,  only  the  data 
involving  the  PA  741  operational  amplifiers  were  analyzed.  The  data  at 
4  MHz  and  15  MHz  are  tabulated  in  Figs.  6.2-1  and  6.2-2,  respectively,  along 
with  the  computed  values  of  1 H^ (— f c » f Q) I • 

The  probability  density  function  associated  with  ) H2 (~f c * f c) l  was  <*e_ 


noted  by  f  (h).  The  procedure  used  to  obtain  approximations  to  f  (h)  was 
n  H 


identical  to  that  followed  for  f  (g) .  Let  Nu  denote  the  number  of  operation- 

G  H 


al  amplifiers  measured.  Also,  let  m^  and  denote  estimates  of  the  mean 


and  k 


th 


central  moment  of  the  random  variable  H.  Finally,  let  6^  and 


denote  estimates  of  the  coefficients  of  skewness  and  kurtosis  for  f„(h). 

H 


If  h,  represents  the  measured  second-order  nonlinear  transfer  function  of 

K 


th 


the  k  operational  amplifier,  then  expressions  for  evaluation  of  m^, 


^H3*  ^H4*  ^H3*  an<*  ^H4  are  identical  i-n  f°rm  to  those  given  by  Eqs.  (6.1-4) 


through  (6.1-9). 


As  was  the  case  for  f  (g) ,  the  x  statistic  was  evaluated  using  the 

G 
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Device  Number 


V2m  (Volts) 


V3m  (Volts>  |H2(-fc,fc)|  (Volts'1) 


1 

.0302 

2 

.0358 

3 

.0238 

4 

.053 

5 

.026 

6 

.0294 

7 

.054 

8 

.0298 

9 

•0337 

10 

.0337 

.1.1 

.0305 

13 

.0365 

14 

.028 

15 

.0332 

16 

.0322 

19 

.0413 

41 

.026 

42 

.034 

43 

.0217 

44 

.024 

45 

.0219 

46 

.023 

47 

.0238 

48 

.0278 

51 

.0265 

52 

.0212 

53 

.0243 

54 

.0224 

55 

.0232 

56 

.0285 

57 

.0278 

58 

.0263 

59 

.0275 

60 

.032 

0.2 

77.5 

0.2 

55.2 

0.2 

124.8 

0.2 

25.2 

0.2 

104.6 

0.2 

81.8 

0.2 

24.2 

0.2 

79.6 

0.2 

62.3 

0.2 

62.3 

0.2 

76.0 

0.2 

53.1 

0.2 

90.2 

0.2 

64.2 

0.2 

68.2 

0.2 

41.5 

0.2 

104.6 

0.2 

61.2 

0.2 

150.2 

0.2 

122.8 

0.2 

147.4 

0.2 

133.7 

0.2 

124.8 

0.2 

91.5 

0.2 

100.7 

0.2 

157.3 

0.2 

119.8 

0.2 

140.9 

0.2 

131.4 

0.2 

87.1 

0.2 

91.5 

0.2 

102.2 

0.2 

93.5 

0.2 

69.1 

Fig.  6.2-2. 


Measurements  and  calculation  of  j H_ (— f  ,f  )| 

2  c  c 

at  15  mHz . 
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relation 
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(pk  -  0.2)2 
0.2 


(6.2-7) 


For  the  data  of  Figs.  6.2-1  and  6.2-2 

Nu  -  34.  (6.2-8) 

H 

2 

Once  again,  the  numerical  value  obtained  for  the  x  statistic  was  compared 
to  a  threshold  value  tabulated  for  a  significance  level  given  by 


Bc=  0.01  (6.2-9) 

and  the  distribution  under  test.  The  cases  for  f  =  4MHz  and  f  =15  MHz 

c  c 

are  discussed  separately. 

Case  1 :  f  =4  MHz 
_ c _ 

Using  the  data  from  Fig.  6.2-1,  the  estimated  mean,  variance,  coefficient 
of  skewness,  and  coefficient  of  kurtosis  were  calculated  to  be 


A  2 
m^  =  2.519  x  10 


=  a 


H2 


H3 


a. 


H4 


-  2.048  x  104 
=  1.480 
=  3.703 


Volts 

Volts 


(6.2-10) 


The  selected  candidates  for  approximating  fu(h)  are  listed  below  along  with 

n 

their  coefficients  of  skewness  and  kurtosis: 

a)  Gamma  (a  =  1),  a ^  =  1.41,  =  6 

b)  Gamma  (a  -  7),  =  0.71,  =  3.75 

c)  Weibull  (8=  2) ,  aj  ■  0.63,  =  3.24 

d)  Log  Normal  (c  =  1/4),  -  0.778,  =  4.096 
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four  distributions. 


Application  of  the  x  test  resulted  in  rejection  of  all 

2 

However,  the  x  test  works  best  when  the  number  of  data  points  is  much 

greater  than  30.  Because  N  =  34,  it  was  decided  to  utilize  the  Kolmogorov  - 

Smirnov  test  which  works  well  when  the  number  of  data  points  equals  30  or 

less  [l6l*  As  before,  the  significance  level  of  the  test  was  set  at  B&=  0.0L 

The  Kolmogorov-Smirnov  test  resulted  in  the  acceptance  of  all  four  candidate 

distributions.  Inserting  the  sample  mean  and  variance  from  Eqs . (6 . 2-10) ,  the 

analytical  expressions  for  the  approximations  to  f  (h)  are  given  by: 

n 


Gamma  (a  =  1) 


(97.65  x  10  6) (h-49.5)exp  {-(9.882  x  10~ 3)  (h-49.5) } ,  h  >  49.5 


fH(h)  * 


Gamma  (a  =  7) 


,  h  <  49.5 

(6.2-11) 


-2 


—  1  8  7 

(4.618  x  10  °)(h+152.9)exp  {-(1.976  x  10  )(h+152.9)},  h<>  -152.9 


fH(h)  ^ 


,  h  <  -152.9 


(6.2-12) 


Weibull  (8  °  2) 


fH(h> 


(4.55  x  10"4)(h-193.1)exp  {-(2.275  x  10~4) (h-193. I)2} , 


0 


» 


h  >  193.1 


h  <  193.1 


(6.2-13) 
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Log  Normal  (C  °  1/4) 


1.596  (h-241.7)"1  exp  {-8 [ £n ( < . 102) (h-241.7))]2},  h  >  241.7 


fH(h) 


,  h  <  41.7. 


(6.2-14) 


The  four  density  functions  are  plotted  in  Figs.  6. 2-3 (a),  (b),(c), 
and  (d) ,  respectively. 

Case  2:  fc_  =  15  M,Hz 

Using  the  data  from  Fig.  6.2-2,  the  estimated  mean,  variance,  coeffi¬ 
cient  of  skewness,  and  coefficient  of  kurtosis  were  calculated  to  be 


*  9.177  x  10  Volts 


Art  A  A  rt 

(°H>  *  MH2  "  1,230  X  10  Volts 


(6.2-15) 


aR3  =  0.048 


a„.  =  2.217. 
H4 


The  selected  candidates  are  listed  below  along  with  their  coefficients  of 
skewness  and  kurtosis: 

a)  Uniform,  ou  =  0,  a  =1.8 

J  4 

b)  Triangular,  =  0,  =  2.4 

c)  Weibull  (8  =  3),  =  0.16,  =  2.73. 

2 

Application  of  the  x  test  resulted  in  acceptance  of  all  three  distribu¬ 
tions.  Inserting  the  sample  mean  and  variance  from  Eqs. (6.2-15) ,  the  analy¬ 
tical  expressions  for  the  three  distributions  are  given  by: 


Fig.  6.2-3.  Approximations  to  f^(h)  at  4  MHz.  (a)  Gamma  (a  =  1_)_» 

(b)  Gamma  (a  =  7),  (c)  Weibull  (6  *  2),  (d)  Log  Normal  ( 


Uniform 


fH(h>  * 


31 .0  <  h  <152.5 


elsewhere 


Triangular 


fH(h) 


(1.36  x  10 


-4 


s  ^  (1 .36  x  10 


-4 


) (h  -  5.86) 
) (177. 7  -h) 


5.86  <  h  <91.8 
91.8  <  h<177.7 


elsewhere 


(6.2-16) 


(6.2-17) 


Weibull  (B  =  3) 


f  00  s 
H 


r. 


039  x  10  ^)(h-81.6)  exp{-(6.8xl0  ^) (h-81. 6) , 


h>81.6 

h  <  81.6 
(6.2-18) 


The  three  density  functions  are  plotted  in  Figs.  6. 2-4 (a),  (b),  and  (c) , 
repectively. 


6.3  Analysis  of  Data  Related  to  Lighting  the  LED 

The  LED  was  observed  to  light  when  the  peak  voltage  of  the  1kHz  tone  at 

V  was  approximately  equal  to  1.5V.  With  f  =  4MHz,  it  was  always  possible  to 
J  c 

light  the  LED.  For  convenience,  the  value  of  the  threshold  at  was  set  at 
A3T  -  vf  =  1.414  Volts.  (4  MHz)  (6.3-1) 

The  rms  reading  at  V^  was  then  equal  to  one  volt  which  was  a  convenient 
level  at  which  to  set  the  threshold. 

At  15  Mhz,  however,  it  was  not  always  possible  to  light  the  LED.  In 
order  to  maximize  the  number  of  available  data  points,  the  value  of  the 
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threshold  at  was  set  at 


A3T  =  0.2  /2  =  0.283  Volts.  (15  MHz) 


(6.3-2) 


The  rms  reading  at  was  then  equal  to  0.2  Volt  which  was  also  a  convenient 
level  at  which  to  set  the  threshold. 

The  amplitude  of  the  carrier  at  V^  which  produced  the  threshold  voltage 
at  V^  was  denoted  by  the  random  variable,  A^.  The  data  obtained  at  4  MHz 
and  15  MHz  are  tabulated  in  ascending  order  in  Figs.  6.3-1  and  6.3-2,  res¬ 
pectively.  The  experiment  involving  lighting  of  the  LED  was  carried  out  at 
the  University  of  Kentucky  while  the  experiment  involving  measurement  of 
the  second-order  nonlinear  transfer  function  was  performed  at  SUNY,  Buffalo. 
The  printed  circuit  boards  at  both  locations  were  designed  to  be  "identical." 
To  verify  whether  this  was  the  case,  several  measurements  were  repeated  at 
SUNY,  Buffalo  using  several  of  the  operational  amplifiers  employed  at  the 
University  of  Kentucky.  At  4 MHz,  the  readings  at  SUNY,  Buffalo  were  approxi¬ 
mately  a  factor  of  1.46  smaller  than  those  at  Kentucky.  At  15  MHz,  the 
readings  at  SUNY,  Buffalo  were  approximately  a  factor  of  1.2  smaller  than 
those  at  Kentucky.  Because  the  second-order  nonlinear  transfer  functions 
were  measured  at  SUNY,  Buffalo,  it  was  decided  to  reduce  the  data  in 
Fig.  6.3-1  by  a  factor  of  1.46  and  the  data  in  Fig.  6.3-2  by  a  factor  of  1.2. 
The  probability  density  function  associated  with  A^  was  denoted  by 

f.  (a.  ) .  The  procedure  used  to  obtain  approximations  to  f.  (a._)  was 
TT  IT  2 

identical  to  that  followed  for  f  (g)  and  f  (h) .  As  before,  the  x  statis- 

G  H 

tic  was  evaluated  using  the  relation. 


2  ^  (p  -  0.2)' 

X  =  na  l  - 


IT  k=l 


0.2 


(6.3-3) 
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Corrected  Data 

A1T  (Volts,  peak)  A^/1.46  A^  (Volts,  nns) 


1.35 

0.92 

1 

1.48 

1.01 

1 

1.50 

1.03 

1 

1.61 

1.10 

1 

1.73 

1.18 

1 

1.82 

1.25 

1 

1.94 

1.33 

1 

2.11 

1.45 

1 

2.16 

1.48 

1 

2.22 

1.52 

1 

2.31 

1.58 

1 

2.32 

1.59 

1 

2.32 

1.59 

1 

2.42 

1.66 

1 

2.56 

1.75 

1 

2.60 

1.78 

1 

2.80 

1.92 

1 

3.22 

2.21 

1 

3.28 

2.25 

1 

3.42 

2.34 

1 

3.76 

2.58 

1 

3.78 

2.59 

1 

3.89 

2.66 

1 

4.19 

2.87 

1 

4.53 

3.10 

1 

5.52 

3.78 

1 

Fig.  6.3-1.  Measurement  of  A.„  at  4WHz 


For  the  data  of 


where  N  denotes  the  number  of  measurements  for  A  . 

IT  1T 

Fig.  6.3-1  and  6.3-2 

N.  =26.  (6.3-4) 

IT 

2 

Once  again,  the  numerical  value  obtained  for  the  x  statistic  was  compared 
to  a  threshold  value  tabulated  for  a  significant  level  given  by 

0  =  0.01  (6.3-5) 

v 

and  the  distribution  under  test.  The  cases  for  f  =  4MHz  and  f  =15  Mhz 

c  c 

are  discussed  separately. 

Case  1 :  f  =4  MHz 
_ c _ 

Using  the  data  from  Fig.  6.3-1  divided  by  the  factor  of  1.46,  the  esti¬ 
mated  mean,  variance,  coefficient  of  skewness,  and  coefficient  of  kurtosis 
were  calculated  to  be 


m  =  1.866  Volts 
A 

(0A)2  =  0A2  =  0.524  (Volts)2 

=  0.899 
A3 

=  3.089. 

A4 


(6.3-6) 


The  selected  candidates  for  approximating  f.  (a1T)  are  listed  below 

IT  1 

along  with  their  coefficients  of  skewness  and  kurtosis: 

a)  Log  Normal  (c  =  1/4),  =  0.778,  =  4.096 

b)  Weibull  (0  =  2) ,  cx3  =  0.63,  a4  =  3.24 

c)  Gamma  (a  =  4),  =  0.89,  =  4.2. 
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Application  of  the  X  test  resulted  in  acceptance  of  all  three 
distributions.  Inserting  the  sample  mean  and  variance  from  Eqs. (6.3-6), 
the  analytical  expressions  for  the  three  distributions  are  given  by: 

Log  Normal  (C  =  1/4) 


-1 


1.596  (a1T  +  3.5)  exp  {-8[£n  ((.192)  (a1T  +  3.5))]  },  alT<>-3.5 


IT 


llvl 


f A  (alT) 
IT 


’  IT 


<-3.5 


(6.3-7) 


Weibull  (B  -  2) 


r 


fA  (alT}  " 
IT 


17.78  (a1T~l  .6)exp{-8.889  (aj  -1.6)  }  ,  _>  1.6 


(6.3-8) 


’  alT  <1,6 


Gamma  (a  =  4) 


^11.71  (a1T  -  .25)4exp{-3.088  (a^  -  .25)},  a1T  >  .25 


f*  (a  =■ 


A  IT 

ait  1 


(6.3-9) 


0 


,  a,™.  <  -25 
IT 


The  three  density  functions  are  plotted  in  Figs.  6.3-3  (a),  (b) ,  and 
(c) ,  respectively. 


Case  2:  f  =  15  MHz 

_ c _ 

Using  the  data  from  Fig.  6.3-2  divided  by  the  factor  of  1.2,  the 
estimated  mean,  variance,  coefficient  of  skewness,  and  coefficient  of 


kurtosis  were  calculated  to  be 


FOUR  MHZ  FREQ 


m  =  1.432  Volts 
A 

(8a)2  -  flA2  ■  1.21  (Volts)2  (6.3-10) 

<SA_  =  1.681 
A3 


v  =  5.931 


The  selected  candidates  for  approximating  f  (a.  )  are  listed  below  along 

A1T  1T 

with  their  coefficients  of  skewness  and  kurtosis: 

a)  Log  Normal  (c  =  1/2),  a3  =  1.75,  =  8.898 

b)  Weibull  (8  =  1) ,  =  2,  =  9 

c)  Gamma  (a  =  1),  *  1.41,  *  6. 

2 

Application  of  the  x  test  resulted  in  acceptance  of  all  three  distributions. 
Inserting  the  sample  mean  and  variance  from  Eqs.  (6.3-10),  the  analytical 
expressions  for  the  three  distributions  are  given  by: 

Log  Normal  (c  =  1/2) 


£a  (aii}  ? 

IT 


.798  (a1T+  l.l)_1exp{-2[iin((.549)(a1T+l.l))]2},  >  -1.1 


(6.3-11) 


’  aiT<  ~lA 


Weibull  (3=1) 


,909  exp{ -(.909) (a  -  .33)}  ,  a1T  1  .33 


(6.3-12) 


Gamma  (a  =  1) 


1.15(alx  +  .12) ,exp{-  (1.071)(a1T+  .12)}  ,  a1T  > 


IT 


fA  (alT)  S' 
A1T  1 


IT 


»  a1T  < 


.12 

(6.3-13) 
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The  three  density  functions  are  plotted  in  Figs.  6.3-4  (a),  (b) , 
and  (c) ,  respectively. 

Figs.  6.3-3  and  6.3-4  illustrate  approximations  to  f  (a  T)  which  were 

IT  1 

obtained  directly  from  the  experimental  data.  Eq. (2.4-11)  was  also  used  to 
predict  the  form  of  fA  (a,_)  based  upon  approximations  to  f  (g)  and 

Aim  li  ^ 

f„(h).  The  cases  for  f  =  4MHz  and  f  =15  MHz  are  discussed  separately. 

H  c  c 


Case  1 :  f  =4  MHz 
c 


At  4  MHz  two  approximations  were  obtained  for  f  (g)  while  four  approxi- 

b 

mations  were  found  for  f  (h) •  The  eight  possible  combinations  are  listed 

H 

below: 


Combination 

fG(g) 

1 

Log  Normal 

(c  =  1/4) 

2 

Log  Normal 

(c  =  1/4) 

3 

Log  Normal 

(c  =  1/4) 

4 

Log  Normal 

(c  =  1/4) 

5 

Gamma  (a  = 

7) 

6 

Gamma  (a  = 

7) 

7 

Gamma  (a  = 

7) 

yto 

Gamma  (  a  =  1) 

Gamma  (a  =  7) 

Weibull  (6  =  2) 

Log  Normal  (c  =  1/4) 
Gamma  (a  =  1 ) 

Gamma  (a  =  7) 

Weibull  (6  =  2) 

Log  Normal  (c  =  1/4) 


8 


Gamma  (a  =  7) 


Each  of  these  combinations  was  used  in  Eq.  (2.4-11)  to  predict  approxima¬ 


tions  to  f  (a  ) .  The  integral  was  evaluated  numerically  and  the  results 
A1T  1 

are  plotted  in  Figs.  6.3-5  and  6.3-6.  The  sharp  peaks  resulted  because  the 

step  size  used  in  the  integration  was  chosen  to  be  too  large. 

2 

A  x  test  was  then  applied  to  each  of  the  eight  approximations  for 

f  (a  )  to  determine  whether  the  corrected  data  from  Fig.  6.3-1  and  the 
A1T  11 

approximations  in  Figs.  6.3-5  and  6.3-6  were  statistically  consistent.  In 

every  case  the  a^  -  axis  was  subdivided  into  5  cells  such  that  the  area 

under  the  approximation  for  each  cell  was  approximately  equal  to  0.2. 

ttl 

Let  the  area  for  the  k —  cell  be  denoted  by  p^.  Also,  let  n^  denote 
the  number  of  corrected  data  which  fell  into  the  k—  cell.  Since  there  are 
26  data  points,  the  relative  frequency  of  occurrence  for  the  k^-  cell  was 
defined  to  be 


2, 


•  •  •  » 


5. 


(6.3-14) 


The  x  statistic  was  then  evaluated  using  the  relation 

2 


2  ?  '  pk> 

X  =  26  l  *  ■ 


(6.3-15) 


k=l 


As  was  done  earlier,  the  numerical  value  obtained  for  the  x  statistic  was 

compared  to  a  threshold  value  tabulated  for  a  significance  level  given  by 

3  =  0.01. 

e 

2 

Application  of  the  x  test  resulted  in  all  eight  approximations  being 


accepted. 


15  MHz 


|c. 


Case  2:  f  = 
c 


At  15  MHz  twc  approximations  were  obtained  for  f  (g)  while  three 

G 

approximations  were  found  for  f  (h)  .  The  six  possible  combinations  are 

rl 

listed  below: 


Combination 

1 

2 

3 

4 

5 

6 


fG(g) 


Gamma  (a  =  3) 
Gamma  (a  =  3) 
Gamma  (a  =  3) 
Gamma  (a  =  4) 
Gamma  (a  =  4) 
Gamma  (a  =  4) 


Vh> 


Uniform 
Triangular 
Weibull  (6  =  3) 
Uniform 
Triangular 
Weibull  (B  =  3) 


Each  of  these  combinations  was  used  in  Eq.  (2.4-11)  to  predict  approximations 
to  f  (a  ) .  The  integral  was  evaluated  numerically  and  the  results  are 

ait  11 

plotted  in  Figs.  6.3-7  and  6.3-8.  As  before,  sharp  peaks  resulted  because 
the  step  size  chosen  for  the  integration  was  too  large. 


A  x  test  was  then  applied  to  each  of  the  six  approximations  for  f^ 

(a^)  to  determine  whether  the  corrected  data  from  Fig.  6.3-2  and  the 

approximations  in  Figs.  6.3-7  and  6.3-8  were  statistically  consistent. 

The  a^-axis  was  again  subdivided  into  5  cells  for  every  case  such  that 

the  area  under  the  approximation  for  each  cell  was  approximately  equal 
2 

to  0.2.  The  x  statistic  was  evaluated  using  Eq. (6.3-15). 

2 

For  a  significance  level  given  by  8^  =  0.01,  the  x  test  resulted 
in  all  six  approximations  being  accepted. 
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Approximations  to  f  (a.  )  at  kHz  obtained  evaluating 

A1T  11 

Eq.  (2.4-11)  with  f G (g)  approximated  by  Gamma  (a=3)  and  f^Ch) 
approximated  by  (a)  Uniform,  (b)  Triangular,  (c)  Weibull  (3=  3) 


s  v 


Fig.  6.3-8.  Approximations  to  f  (a.  )  at  *5  MHz  obtained  by  evaluating 

IT  1 

Eq.  (2.4-11)  with  f  (g)  approximated  by  Gamma  (a=4)  and  f  (h) 

O  n 

approximated  by  (a)  Uniform, (b)Triangular ,  (c)  Weibull  (8=3 


^  «/*  «."*  k'*»  >,**■  *S  -  *  V‘  -  ‘  ■"*  •'*  *'• 


VII.  Conclusions 


The  LED  experiment  served  to  demonstrate  the  following  results: 

1)  Crosstalk  in  random  cable  bundles  can  lead  to  serious  degradation 
of  system  performance. 

2)  Weak  nonlinearities  which  exist  in  the  amplification  region  of 
operational  amplifiers  are  capable  of  generating  significant  in-band  inter¬ 
ference  by  demodulating  out-of-band  amplitude  modulated  signals. 

3)  The  linear  transfer  function  which  characterizes  crosstalk  in 
random  cable  bundles,  the  second-order  nonlinear  transfer  function  which 
characterizes  the  weak  nonlinearities  in  the  amplification  region  of  opera¬ 
tional  amplifiers,  and  the  susceptibility  level  of  the  three-stage  amplifier 
all  display  random  behavior. 

4)  Well-known  probability  density  functions  can  be  used  to  approximate 
the  crosstalk  linear  transfer  function,  the  second-order  nonlinear  transfer 
function  of  operational  amplifiers,  and  the  susceptibility  level  of  the 
three-stage  amplifier. 

5)  Given  probability  density  functions  for  the  crosstalk  linear  trans¬ 
fer  function  and  the  second-order  nonlinear  transfer  function  of  operational 
amplifiers,  it  is  possible  to  predict  the  probability  density  function  for  the 
susceptibility  level  of  the  three-stage  amplifier. 

6)  The  predicted  probability  density  function  for  the  susceptibility 
level  of  the  three-stage  amplifier  can  be  approximated  by  well-known  proba¬ 
bility  density  functions  that  are  statistically  consistent  with  experimental 
data.  Although  not  included  in  this  report,  these  approximations  are  useful 
for  determining  the  probability  of  EMI  and/or  compatibility. 
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1 4  provided  to  ESP  Program  Office*  (PPa)  and  othei t  ESP 
element*.  The  principal  technical  mi**lon  an. ea*  one 
communication*,  electromagnetic  guidance  and  control,  sur¬ 
veillance  of  ground  and  aerospace  object*.  Intelligence  data 
collection  and  handling,  Information  system  technology. 
Ionospheric  propagation,  solid  state  sciences,  microwave 
physic*  and  electronic  reliability,  maintainability  and 
compatibility. 
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